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The  potential  for  optical  communication  through  low-visibility  scattering 
atmospheres  is  the  underlying  concern  of  this  investigation.  To  access  that 
potential,  the  effects  of  the  channel  upon  the  optical  signal  propagating  through 
it  must  be  determined.  The  work  reported  here  is  primarily  concerned  with 
one  of  those  effects:  the  time  dispersion,  or  multipath  spread,  of  the  trans- 
mitted signal. 

The  investigation  was  primarily  an  experimental  one,  utilizing  the  existing 
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between  the  MIT  campus  and  a field  station  located  at  an  Air  Force  site  in 
suburban  Boston.  Experiments  were  performed  at  visible  and  near  infrared 
wavelengths,  with  highly  collimated  transmitting  beams.  Measurements  of 
the  multipath  spread  for  various  receiver  fields  of  view  and  at  various  off- 
axis  angles  were  made. 

Little  multipath  and  angular  spread  was  observed  within  the  realm  of 
atmospheric  conditions  for  which  a detectable  signal  could  be  obtained.  For 
optical  thicknesses  less  than  ten,  the  on-axis  measurements  are  consistent 
with  the  hypothesis  that  the  received  signal  is  dominated  by  unscattered  radia- 
tion; however,  the  off-axis  measurements  suggest  that  the  scattered  field 
itself  may  be  only  slightly  spread  in  time  and  angle. 

Guided  by  the  observed  narrowness  of  the  angular  spectrum,  an  approxi- 
mation to  the  linear  transport  equation  was  developed.  The  resulting  equation 
is  more  amenable  to  solution  than  is  the  full  transport  equation.  Here  it  is 
solved  for  isotropic  scatter.  For  such  scattering  it  is  found  that  substantial 
spreading  in  both  time  and  angle  will  occur  when  the  optical  thickness  becomes 
appreciably  greater  than  one.  This,  in  combination  with  the  data  collected, 
suggests  that  nonisotropic  scattering  encountered  in  the  atmosphere  is 
responsible  for  the  absence  of  spreading.-^ 
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EVALUATION 


1.  This  report  is  the  Final  Report  on  the  contract.  It  covers 
research  done  on  optical  communication  channels  during  the  15 
month  period  1 July  1975  to  30  September  1976.  The  objective 
of  the  research  is  to  determine  the  effect  of  low-visibility 
scattering  atmospheres  upon  optical  signals.  There  exists  the 
potential  for  improving  the  operability  of  1 i ne-of-s i ght  optical 
communication  links  by  exploiting  the  scattered  component  of  the 
received  optical  signal  in  addition  to  the  unscattered  part. 
Experiments  were  conducted  under  this  contract  to  determine  the 
basic  properties  of  the  scattered  part  of  the  transmitted  optical 
si gnal  . 

2.  The  above  work  is  of  value  since  it  provides  basic  knowledge 
on  optical  communication  channels  under  low-visibility  conditions 
which  will  allow  the  scattered  part  of  the  transmitted  signal  to 
be  usefjl  in  communication  links  in  addition  to  the  unscattered 
part.  This  research  may  make  the  development  of  new  more  sensitive 
optical  receivers  possible  which  will  increase  the  range  of  optical 

communication  channels  designed  for  USAF  use. 
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SUMMARY 


The  potential  for  optical  communication  through  low-visibility  scattering 
atmospheres  is  the  underlying  concern  of  this  investigation.  To  access  that 
potential,  the  effects  of  the  channel  upon  the  optical  signal  propagating  through 
it  must  be  determined.  The  work  reported  here  is  primarily  concerned  with 
one  of  those  effects:  the  time  dispersion,  or  multipath  spread,  of  the  trans- 

mitted signal. 

The  investigation  was  primarily  an  experimental  one,  utilizing  the  existing 
propagation  facility  operating  over  a thirteen  kilometer  experimental  path 
between  the  MIT  campus  and  a field  station  located  at  an  Air  Force  site  in 
suburban  Boston.  Experiments  were  performed  at  visible  and  near  infrared 
wavelengths,  with  highly  collimated  transmitting  beams.  Measurements  of 
the  multipath  spread  for  various  receiver  fields  of  view  and  at  various  off- 
axis  angles  were  made. 

Little  multipath  and  angular  spread  was  observed  within  the  realm  of 
atmospheric  conditions  for  which  a detectable  signal  could  be  obtained.  For 
optical  thicknesses  less  than  ten,  the  on-axis  measurements  are  consistent 
with  the  hypothesis  that  the  received  signal  is  dominated  by  unscattered  radia- 
tion; however,  the  off-axis  measurements  suggest  that  the  scattered  field 
itself  may  be  only  slightly  spread  in  time  and  angle. 

Guided  by  the  observed  narrowness  of  the  angular  spectrum,  an  approxi- 
mation to  the  linear  transport  equation  was  developed.  The  resulting  equation 
is  more  amenable  to  solution  than  is  the  full  transport  equation.  Here  it  is 
solved  for  isotropic  scatter.  For  such  scattering  it  is  found  that  substantial 
spreading  in  both  time  and  angle  will  occur  when  the  optical  thickness  becomes 
appreciably  greater  than  one.  This,  in  combination  with  the  data  collected, 
suggests  that  nonisotropic  scattering  encountered  in  the  atmosphere  is 
responsible  for  the  absence  of  spreading. 
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CHAPTER  I.  INTRODUC  TION 


The  apparent  limitations  imposed  upon  atmospheric  optical  communication 
systems  by  low-visibility  weather  conditions  (haze,  rain,  fog,  snow,  etc.)  have 
been  widely  noted.  In  particular,  for  a conventional  line-of-sight  link,  whose 
receiver  responds  only  to  the  unscattered  light  arriving  along  the  boresight 
from  the  transmitter,  the  major  effect  of  low-visibility  atmospheric  propaga- 
tion is  to  increase  the  path  loss  by  a factor  of  exp(-cL),  where  L is  the  path 
length  and  a is  the  extinction  coefficient  of  the  atmosphere  at  the  wavelength 
of  interest.  Although  the  relationship  between  a and  visibility  is  tenuous,  at 
best,  some  appreciation  for  the  magnitude  of  the  attenuation  is  provided  by  the 

commonly  accepted  approximation  that,  at  visible  wavelengths,  exp(-cL)  will 

_ •> 

be  on  the  order  of  10  ^ when  the  operating  range  1,  equals  the  visibility  [l,  Z\. 

At  most  visible  wavelengths,  and  at  many  infrared  wavelengths,  the  extinc- 
tion coefficient  c is  due  in  large  measure  to  scattering  rather  than  absorption 
[3,4].  Thus,  there  exists  the  potential  for  improving  the  operability  of  line-of- 
sight  optical  links  by  exploiting  the  scattered  component  of  the  received  optical 
field,  in  addition  to  the  unscattered  component  [5-8],  The  extent  to  which  utili- 
zation of  scattered  light  will  significantly  improve  low- visibility  communication 
performance  depends  upon  the  energy  contained  in  the  scattered  component  of 
the  received  field  and  upon  the  way  in  which  this  energy  is  divided  between  the 
statistical  degrees  of  freedom  of  the  field.  These  quantities  are,  in  turn,  re- 
lated to  the  characteristics  of  the  scattering  particles  and  to  the  geometry  of 
the  channel  in  a very  complicated  way.  Nevertheless,  some  basic  features  can 
be  discerned. 

It  is  reasonable  to  assume  that  when  multiple-scattering  effects  predomi- 
nate, as  it  is  likely  they  will  in  a low-visibility  atmospheric  channel,  the  opti- 
cal complex  field  envelopes  received  on  two  orthogonal  polarizations  will  be 
independent  identically-distributed  zero-mean  complex-Gaussian  random  pro- 
cesses whose  statistics  are  completely  characterized  by  a space-time  correla- 
tion function  [9-11  J.  When  the  particle  scattering  function  is  sharply  peaked 
in  the  foward  direction,  various  analytic  techniques  (which  have  been  used  to 
study  wave  propagation  through  atmospheric  turbulence)  may  be  applied  to  ex  - 
plicitly  determine  the  desired  space-time  correlation  function  [12,  1 3 J.  At 
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present,  however,  multiple-scattering  theory  is  too  complicated  to  permit 
such  an  oplicit.  determination  in  the  general  case,  wherein  the  forward- 
scattering assumption  is  no  longer  applicable  [9-11  J. 

Fortunately,  previous  experience  with  radio-frequency  scattering  chan- 
nels suggests  that  detailed  knowledge  of  the  channel  correlation  function  may 
not  be  necessary.  In  particular,  specification  of  "physical"  parameters  such 
as  multipath  spread  (time  dispersion),  Doppler  spread  (frequency  dispersion), 
and  angular  spread  (spatial-frequency  dispersion)  can  suffice  for  both  the 
design  of  receiver  structures  and  performance  calculations  [5,  14 J.  Indeed, 
based  on  reasonable,  but  as  yet  unverified,  statistical  assumptions,  two  pos- 
sible receiver  structures  for  exploiting  the  information  contained  in  the  scat- 
tered field  have  been  identified  [5,  6,  8 J. 

The  first  structure  is  a wide  field-of-view  direct  detection  system  that 
counts  photon  arrivals  in  the  entire  time -frequency- spatial  frequency  region 
over  which  scattered  signal  light  is  expected.  The  performance  obtainable 
with  this  receiver  will  depend  critically  upon  the  development  of  improved 
wide-angle  narrow-band  apical  filters  to  discriminate  against  background 
light,  or  operation  at  wavelengths  where  such  light  is  negligible.  Under  favor- 
able conditions,  the  direct-detection  performance  can  approximate  that  of  the 
optimum  one-shot  quantum -mechanical  receiver.  The  second  candidate 
receiver  is  a spatial  phase-compensation  (mode-compression)  receiver,  sim- 
ilar to  systems  previously  recommended  for  atmospheric  turbulence  (clear- 
weather)  applications  |l5-17  J,  but  with  many  more  phase-controllable  elements. 
Ihe  theoretical  utility  of  mode  compression  rests  on  the  assumptions  that  chan- 
nel fluctuations  occur  slowly  relative  to  the  desired  information  rate,  and  that 
the  uncompensated  system  is  background-limited.  The  practicality  of  the 
phase-compensation  receiver  will  depend  upon  the  ease  with  which  a large  num- 
ber of  phase  elements  may  be  reliably  estimated  and  controlled. 

In  the  context  of  the  foregoing  receiver  structures,  the  importance  of  the 
channel  spread  parameters  (multipath,  Doppler,  etc.)  is  apparent,  and  the  need 
for  experimental  measurements  manifest.  The  degree  to  which  these  receivers 
will  extend  the  operability  of  atmospheric  links  depends  overwhelmingly  on  the 
degree  to  which,  under  low-visibility  conditions,  the  energy  contained  in  the 
scattered  component  of  the  received  field  substantially  exceeds  the  energy  con- 
tained in  the  ur  scattered  component.  For  example,  it  can  be  argued  that  for  a 
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beam  source  embedded  in  an  infinite  lossless  homogeneous  scattering  medium, 
the  received  energy  at  path  length  L will,  as  visibility  decreases,  ultimately 
reduce  to  tne  level  of  a free -space  omnidirectional  radiator.  Whether  or  not 
the  scattered-field  energy  still  dominates  the  unscattered  energy  by  a wide 
margin  in  the  more  typical  case  of  a partially  absorbing  inhomogeneous  scat- 
tering medium  with  boundaries,  remains  to  be  determined  (cf.  Chap.  II). 

Given  '.hat  there  is  sufficient  energy  in  the  scattered  field  to  warrant  the 
use  of  a low-visibility  (wide  field-of-view  or  mode-compression)  receiver,  the 
ease  with  which  such  a system  may  be  implemented  depends,  primarily,  on  the 
angular  spread  of  the  scattered  light.  This  angular  spread  must  be  accommo- 
dated by  the  optical  filter  in  the  wide  ficld-of-view  receiver.  This  angular 
spread,  measured  in  units  of  diffraction-limited  fields  of  view,  determines  the 
required  number  of  phase  controllers  in  the  mode-compression  receiver.  In 
either  case,  the  implementation  burden  is  an  increasing  function  of  angular 
spread. 

Finally,  assuming  that  a substantial  amount  of  scattered  energy  can  be 
collected  over  an  acceptable  (vis-a-vis  implementation)  field  of  view,  the  pri- 
mary limitation  on  the  two  candidate  low-visibility  receivers  may  be  due  to 
multipath  spread.  In  particular,  to  achieve  data  rates  in  excess  of  the  recip- 
rocal multipath  spread  for  the  full  (uncompensated)  receiver  field  of  view  requires 
an  additional  level  of  system  complexity  to  overcome  intersymbol  interference. 
Moreover,  because  multipath  spreads  in  excess  of  one  microsecond  have  been 
measured  in  a cloud  propagation  experiment  1 1 8 , 1 9 J,  the  data  rate  constraint 
on  simple  receivers  may  be  severe. 

The  preceding  considerations  have  led  us  to  make  the  nature  of  the  scattered- 
energy  distribution  a central  concern  of  our  ongoing  experimental  and  theoretical 
research  on  the  possibilities  for  improved  optical  communication  in  low- 
visibility  environments  [ 6,  20  J.  The  objectives  of  the  experimental  program 
have  been  to: 

1)  measure  the  ratio  of  scattered  to  unscattered  radiation  collected  by  a 
receiver  as  a function  of  its  field  of  view,  and  meteorological  conditions; 

2)  measure  the  angular  spectra  of  pulsed  and  cw  optical  signals  as  functions 
of  meteorological  conditions; 

3)  measure  Ihe  multipath  spread  of  pulsed  optical  signals  as  a function  of 
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meteorological  conditions;  and 

-U  measure  tin-  Doppler  spread  of  cw  optical  signals  as  a function  of 
meteorological  conditions. 

The  approach  taken  in  the  theoretical  work  has  been  to  solve  the  linear- 
transport  equation,  subject  to  approximations  motivated  by  experimental 
results,  for  the  parametric  dependence  of  angular  spread,  beam  spread,  and 
multipath  spread  on  the  total  extinction,  albedo,  particle  scattering  function, 
and  source  geometry. 

We  shall  report  here  on  the  experimental  and  theoretical  studies  of  multi- 
path  spread.  These  studies  were  funded  by  the  Air  Force  Electronic  Systems 
Division  under  Contract  F 19628-76-C-0054  1 2 1 I. 


4 
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CHAPTER  II.  EXPERIMENTS 


As  previously  noted,  multiple-scattering  theory  (in  its  present  state)  does 
not  provide  an  explicit  determination  of  the  critical  spread  parameters  which 
characterize  the  low-visibility  atmospheric  channel.  We  have  therefore  been 
pursuing,  under  National  Science  Foundation  and  Air  Force  Electronic  Systems 
Division  funding  [21, 22],  an  experimental  program  aimed  at  measuring  these 
spread  parameters.  In  this  chapter,  our  multipath  measurements  will  be  dis- 
cussed beginning  with  a description  of  the  propagation  geometry  and  measure- 
ment apparatus. 


A.  Propagation -Link  Specifications 

Tlie  multipath  data  reported  in  the  sequel  was  obtained  on  a propagation 
facility  established  under  NSF  funding  over  a 13.  (>  km  line-of-sight  path  between 
the  optical  communications  laboratory  atop  the  Earth  Sciences  building  at  M.I.T. 
and  the  Air  Force  Millimeter- Wave  Observatory  atop  Prospect  Hill  in  Waltham 
Mass.  The  latter  location  has  been  made  available  to  us  through  the  generous 
cooperation  of  Dr.  E.  E.  Altshuler  and  his  group  at  the  Deputy  for  Electronic 
Technology.  A profile  of  the  terrain  along  the  path  is  shown  in  Fig.  1,  and 
photographs  of  the  views  along  the  path  from  eacli  terminal  are  shown  in  Figs.  2 
and  3. 

This  propagation  facility  comprise  ./four  links  utilizing  continuous-wave 
helium  neon,  pulsed  ruby,  and  pulsed  holmium  lasers.  Multipath  measurements 
were  made  exclusively  with  the  ruby  link.  However,  important  corroborating 

evidence  was  obtained  from  holmium-link  angular  spectra.  These  two  links 

\ 

are  described  below. 

1 . Ruby  Link 

A block  diagram  of  the  ruby  link  is  shown  in  Fig.  4.  The  transmitter  is 
a Q-switched  Spacerays  ruby'-laser,  which  is  located  at  M.  I.  T.  This  laser 
delivers  pulses  of  nominal  0.  5 J energy  and  70  nsec  duration  at  a rate  of  10  per 
minute.  The  transmitter  beam  divergence  is  controlled  by  a beam-expansion 
telescope.  Initial  measurements  wore  made  with  a 1.  8 mRad  full  width  to  half 
maximum  (FWHM)  beam  divergence;  the  telescope  was  later  modified  to  yield 
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Fig.  1.  Profile  of  terrain  along  propagation  path. 
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A photograph  of  the  ruby-link  receiver  is  shown  in  Fig.  5.  The  optical 
subsystem  consists  of  a modified  Zoomar  telescope  (127  cm"  effective  aper- 
ture, f/3.  5,  Schmidt-Cassegrain)  mounted  on  a pair  of  rotary  tables  for  azi- 
muth and  elevation  control.  These  tables  are  equipped  with  electronic  resolvers 
which  permit  absolute  aiming  and  resetting  to  0.  1 mRad.  Pinholes  may  be 
inserted  into  the  focal  plane  to  select  fields  of  view  from  0.  1 mRad  to  20  mRad 
full  cone  angle.  Kodak  Wratten  neutral  density  filters  are  employed  (in  clear- 
to  modest-visibility  weather)  at  this  point  in  the  system  to  prevent  overloading 
the  photodetector. 

After  the  focal  plane,  the  ruby  light  is  recollimated  and  transmitted  through 

a dichroic  beamsplitter  (which  enables  the  telescope  optics  to  be  shared  by  the 

holmium-link  receiver)  and  a lnm  linewidth  694.  3 nm  wavelength  interference 

filter  onto  the  photodetector.  The  ruby  receiver  employs  an  EMI  9558B  S-20 

cathode  photomultiplier  which  has  a responsivity  of  0.  018  A/W  at  694.  3 nm 

and  a current  gain  of  10^  (as  operated).  The  detector  has  a mean  anode  dark 
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current  (at  room  temperature)  of  10  A,  and,  including  the  effects  of  the  post- 
detection electronics,  yields  single-photon  events  with  a 30  nsec  rise-fall  time. 

The  post-detection  signal  processing  for  the  ruby  link  consists  of  a Tek- 
tronix Mode)  466  100  MHz  storage  oscilloscope,  and  a P.  A.  R.  Model  162 
boxcar  integrator  with  Model  9002  X-Y  recorder.  Because  of  the  extremely 
low  duty  cycle  of  the  ruby  laser,  the  normal  data-collection  procedure  is  to 
view  each  individual  pulse  on  the  storage  oscilloscope  and  record,  by  hand,  the 
peak  voltage  and  time  duration.  Permanent  photographic  records  are  made  of 
unusual  pulse  shapes.  Under  severe  weather  conditions,  for  which  oscilloscope 
observations  are  insufficiently  sensitive,  boxcar  averaging  of  5-20  pulses  is 
performed,  and  the  results  are  plotted  by  the  X-Y  recorder.  Timing  informa- 
tion for  both  the  oscilloscope  and  the  boxcar  integrator  is  provided  by  a 450  MHz 
radio  link,  which  is  described  in  section  A.  3. 

2.  Holmium  Link 

The  holmium  link,  shown  in  block  diagram  in  Fig.  6,  is  built  around  a 
2.06  p wavelength  Ho:YLF  pulsed  laser  developed  by  the  Center  for  Materials 
Science  and  Engineering  at  M.  I.  T.  for  use  in  atmospheric  propagation  experi- 
ments [2 3 J.  The  2.06  p wavelength  has  several  virtues  for  optical  communication 
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Fig.  6.  Holmium-link  block  diagram. 


in  the  atmosphere.  Mast  notably,  because  this  wavelength  is  strongly  absorbed 
by  the  cornea  and  aqueous  humor  of  the  eye,  it  is  relatively  "eyesafe"  in  that 
the  energy  brought  to  focus  on  the  retina  is  reduced  considerably  compared,  for 
example,  to  the  1.06  p wavelength  radiation  of  a NckYAG  laser.  The  holmium 
laser  used  in  the  propagation  studies  is  located  at  M.  I.  T.  , and  operates  in  the 
long-pulse  mode.  It  delivers  pulses  of  nominal  energy  0.  1 J and  150  psec  dura- 
tion at  a rate  of  5 per  second.  Early  measurements  were  made  using  the 
2 mRad  FWHM  natural  beam  divergence  of  the  laser.  Later  measurements 
were  made,  using  a diverging  lens,  with  a 5.  5 mRad  FWHM  beam  divergence. 

The  holmium-link  receiver  shares  its  principal  optical  components  and 
post-detection  signal  processing  with  the  ruby-link  receiver.  Specifically,  the 
holmium  light  is  brought  to  focus  by  the  Zoomar  optics,  with  the  same  azimuth, 
elevation,  and  field  of  view  control  cited  in  section  A.  1.  In  clear  to  modest 
visibility  weather,  Kodak  Wratten  neutral  density  filters  (which  have  been  cali- 
brated at  2.  0 p wavelength)  are  inserted  in  the  focal  plane  to  prevent  overloading 
the  photodetector.  After  the  focal  plane,  the  holmium  light  is  recollimated, 
and  reflected  by  the  dichroic  beamsplitter  through  a 33.  5 nm  linewidth,  2.06  p 
wavelength  interference  filter  off  a 6:  1 demagnification  elliptical  mirror  and 
onto  the  photodetector.  The  holmium  receiver  uses  a Judson  J-12  InAs  photo- 
diode mounted  in  a liquid  nitrogen  dewar.  This  detector  has  a responsivity  of 
0.12  A/W  at  2.  06  p and  a 1 . 0 psec  time  constant.  Its  output  is  coupled  to  a low- 
noise  preamplifier  whose  nominal  current  gain  is  10*\  and  whose  root  mean 
square  (RMS)  feise  current  is  1 mA  (both  measurements  made  with  a 50  ohm  load). 
The  preamplifier  output  may  be  monitored  on  the  Tektronix  466  oscilloscope. 
Normal  data  taking  procedure,  however,  involves  boxcar  averaging  of  20  pulses 
and  X-Y  recording;  longer  pulse  trains  are  averaged  in  the  poorest  visibility 
weather.  Synchronization  is  provided  by  the  450  MHz  radio  link  described 
below. 

3.  Synchronization  Link 

Timing  information  for  both  the  ruby  and  holmium  links  is  provided  by  a 
450  MHz  radio  link  [24  J,  which  consists  of  two  Motorola  U44  transceivers  mod- 
ified for  pulsed  operation.  Trigger  pulses  from  either  laser  are  transmitted 
from  Cambridge  to  Waltham.  A transmitter  output  power  of  approximately  one 
watt  and  six-element  Yagi  antennas  at  both  ends  of  the  line-of-sight  path  give  a 
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signal-to-noise  ratio  in  excess  of  20  dll  prior  to  pulse  regeneration.  ITie  use 
of  regulated  power  supplies  at  l>oth  t>  rminals  renders  long-term  timing  drift 
negligible  after  a 30  minute  warmup  period.  The  short-term  timing  jitter  has 
been  measured  at  ±100  nsec  or  less  (see  Fig.  7).  For  the  ruby  link,  the  sync 
pulse  arrives  approximately  3 psec  before  the  optical  pulse;  for  the  holmium 
link,  the  sync  pulse  arrives  approximately  240  psec  before  the  optical  pulse. 

B.  Propagation-Link  Predictions 

The  impetus  for  propagation  experiments  derives  from  the  dearth  of  avail- 
able theoretical  or  experimental  results  for  the  crucial  atmospheric  spread 
parameters.  Perhaps  the  only  prior  effort  at  characterizing  a low-visibility 
optical  communication  channel  was  the  work  done  on  cloud  propagation  chan- 
nels by  Bucher  et  al.  at  the  M.  I.  T.  Lincoln  Laboratory  [ 1 8,  19,25].  This 
program,  which  sought  to  determine  the  viability  of  satellite  optical-frequency 
downlinks  in  the  presence  of  cloud  cover  at  the  ground  station,  included  a cloud- 
propagation  experiment  and  a Monte-Carlo  cloud-propagation  computer  simu- 
lation. In  this  section  we  shall  apply  the  scaling  laws  obtained  in  the  foregoing 
computer  simulation  in  an  attempt  to  characterize,  a priori,  the  behavior  of 
our  propagation  link  under  low-visibility  conditions.  We  recognize  that, 
because  conditions  present  on  our  link  violate  a number  of  premises  made  in 
the  Lincoln  Laboratory  simulation,  there  is  no  guarantee  that  the  following  pre- 
dictions will  be  borne  out  in  the  experiment.  We  feel,  however,  that  the 
Lincoln  Laboratory  work  provides  a useful  setting  in  which  to  discuss  our 
experimental  results. 

1.  Lincoln  Laboratory  Simulation 

We  may  succinctly  summarize  Bucher's  Monte  Carlo  cloud  simulation  as 
follows.  The  atmosphere  was  modeled  as  a plane-parallel  slab  of  homoge- 
neously distributed  scatterers.  The  slab  had  thickness  L and  was  infinitely 
wide.  The  particles  were  assumed  to  be  nonabsorbing  (unity  albedo),  with  a 
common  scattering  function  which  was  taken  variously  to  be  a Mie  pattern,  a 
Henyey-Greenstein  function,  or  a "megaphone"  (unique  scattering  angle).  The 
parameters  of  each  of  these  scattering  functions  were  adjusted  to  yield  an  aver- 
age cosine  of  the  scattering  angle,  (cos  6),  of  0.  875. 
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For  a ray  source  normally  incident  on  the  atmospheric  slab,  the  distribu- 
tion of  transmitted  rays  emerging  from  the  slab  was  determined  by  a Monte 
Carlo  procedure.  Significantly,  it  was  found  that  the  output  distribution  in 
space,  angle,  and  time  was,  for  thick  enough  slabs,  independent  of  which  scat- 
tering function  was  employed  in  the  simulation.  Bucher  argues  that  this  behav- 
ior is  characteristic  of  propagation  when  multiple-scattering  predominates. 

He  deduces  the  following  scaling  laws  for  the  multiple-scattering  limit: 

i)  the  fraction  of  the  incident  radiation  that  is  transmitted  through  the 
slab,  T,  satisfies 

T~  1.69/(rd+1.42),  (1) 

where  = cL(l-(cos  G)); 

ii)  the  light  transmitted  through  the  slab  exits  as  a circular  spot  of  medi- 
an diameter,  d , which  satisfies 

ds  = 1.56  L/t®'  °7;  (2) 

iii)  the  angular  spectrum  of  the  transmitted  radiation  is  quite  broad,  it 
is  concentrated  within  a tt/4  Rad  full  cone-anph  field  of  view  with  no 
sharp  peaks  to  the  distribution; 

ivl  the  RMS  time  dispersion  of  the  transmitted  radiation,  o(,  satisfies 
trt  ~ 0.  64  L(rd)°- 82/c,  (3) 

where  c is  the  speed  of  light. 

The  foregoing  results  apply  to  slabs  for  which  s 3.  Bucher  attaches  the 
term  diffusion  regime  to  propagation  conditions  which  satisfy  the  above  crite- 
rion. His  justification  for  this  terminology  is  that  computer  simulation  of  pure 
diffusion  (isotropic  scatter)  in  three  dimensions  yielded  scaling  laws  for  dR  and 
that  were  similar  to  Eqs.  (2,  3)  when  compared  as  functions  of  (rather  than 
t).  The  parameter  is  therefore  the  diffusion  thickness  of  the  slab,  i.  e.  , it 
is  the  usual  optical  thickness,  t = oL,  modified  by  the  predilection  of  the  par- 
ticles to  scatter  in  the  forward  direction. 
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2.  Transmission- Fraction  Analysis 

Let  us  apply  Bucher's  scaling  laws  to  the  Cambridge-to- Waltham  propaga- 
tion path.  We  shall  suppose  that  these  laws  are  valid  for  all  values  of  (cos  0). 
Moreover,  we  shall  not  attempt  to  account  for  atmospheric  inhomogeniety, 
particle  absorption,  or  boundary  (ground)  losses,  that  are  obviously  present 
to  some  extent  in  our  experiment,  but  were  not  induced  in  the  computer  simu- 
lation. Our  primary  concern  herein  will  be  10  delineate  two  asymptotic  propa- 
gation regimes  for  this  path:  the  clear-weather  regime  wherein  unscattered 

radiation  dominates  the  received  signal,  and  the  very  low  visibility  regime 
wherein  diffuse  radiation  (governed  by  the  scaling  laws  of  section  B.  1)  domi- 
nates the  received  signal.  In  particular,  we  shall  determine  the  crossover 
optical  thickness  at  which  these  two  asymptotes  are  equal.  However,  because 
of  factors  neglected  in  this  analysis,  we  cannot  argue  that  the  calculated  cross- 
over optical  thickness  (which  under  some  circumstances  does  not  exist)  accu- 
rately characterizes  the  Cambridge  to  Waltham  path.  Indeed,  we  cannot  deny 
the  possibility  that  there  is  another  mode  of  scatter  propagation  which  renders 
the  above  crossover-thickness  calculation  valueless.  Putting  these  caveats 
aside,  let  us  begin  by  considering  the  unscattered  radiation. 

The  fraction,  T , of  the  transmitter  signal  that  reaches  the  plane  contain- 
ing the  receiver's  entrance  pupil  as  unscattered  light  is  given  by  the  Bouguer- 
Lambert  extinction  law 

Tu  = exp(-T).  (4) 

Of  the  unscattered  light  reaching  the  receiver  entrance -pupil  plane,  only  a 
fraction  is  intercepted  by  the  entrance  pupil,  where  A^  is  the  effective 

pupil  area  and  0t  is  transmitter  beam  divergence  (full  cone  angle).  Thus  Tr, 
the  fraction  of  the  transmitter  signal  that  is  intercepted  as  unscattered  radiation 
by  the  receiver  pupil,  satisfies 

= Ar  exp(-T)/(0tL)2.  (5) 

When  the  transmitter  and  receiver  are  in  boresight  alignment,  essentially  all 
of  the  unscattered  radiation  entering  the  receiver  pupil  will  pass  through  the 
focal-plane  pinhole  and  onto  the  photodetector  if  the  receiver  field  of  view 
(determined  by  the  pinhole)  exceeds  the  nominal  clear-weather  blur  circle  of 
the  telescope.  This  blur  circle  would  have  a full  cone  angle  of  the  order 
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of  X/Ar  “ Rad  at  wavelength  X for  a diffraction-limited  instrument.  Tele- 
scope aberrations,  atmospheric  turbulence,  and  platform  vibration  may  typi- 
cally increase  the  blur  circle  full  cone  angle  to  0.  1 mRad. 

l et  us  now  consider  the  diffusion  regime.  When  the  visibility  is  low 
enough  to  satisfy  t(J  > 3,  the  fraction,  T,  of  the  transmitter  signal  (hat  reaches 
the  plane  containing  the  receiver's  entrance  pupil  will  be  assumed  to  satisfy 
Eq.  (1).  In  Fig.  8 we  have  plotted  Tu  and  T versus  optical  thickness,  t = aL, 
for  various  values  of  (cos  9).  The  T curves  all  begin  at  Td  = 3.  Figure  8 
shows,  in  essence,  what  was  suggested  in  Chap.  I,  i.  e.  , for  a lossless  scat- 
tering medium  the  transmission  of  scattered  radiation  may  enormously  exceed 
the  unscattered  transmission. 

Of  the  scattered  light  reaching  the  receiver  entrance-pupil  plane,  the 
fraction  that  is  collected  by  the  pupil  will  be  assumed  to  be  A /d2,  where  d 
satisfies  Eq.  (2).  Thus,  the  fraction  of  the  transmitter  light  thaUin  the  diffusion 
regime)  is  collected  by  the  receiver  pupil,  Tr,  is 

Tr*  0.696  Ar(Td)°-14/L2(Td  + 1.42).  (6) 

In  b ig.  9 we  have  plotted  T^  and  Tr  versus  optical  thickness  for  various  values 
of  'cos6)  using  Ar  = 127  cm2,  6t  = 10  mRad,  and  L = 13.6  km.  Evidently, 
for  sufficiently  large  (cose)  values  (e.  g.  , (cos  0>  = 0.  8 in  Fig.  9)  there  will 
not  be  a crossover  optical  thickness  at  which  the  unscattered-light  and  the 
diffusion -regime  receiver-pupil  transmission  fractions  are  equal,  i.  e.  , for 
sufficiently  large  < cos  0)  we  have  T^  < T*  at  T(1  = 3.  However,  when  there  is 
a crossover  optical  thickness  in  Fig.  9,  it  occurs  in  the  vicinity  of  t = 1 1 , an 
optical  thickness  for  which  the  diffusion-regime  receiver-plane  transmission 
fraction  greatly  exceeds  the  corresponding  transmission  fraction  for  unscattered 
light,  c.f.  big.  8.  This  is  due  to  the  enormous  spatial  spreading  ot  the 
multiply-scattered  light  as  opposed  to  the  small  spatial  extent  of  the  highly 
collimated  unscattered  beam.  The  crossover  point  shifts  to  higher  optical 
thicknesses  (lower  visibilities)  when  the  effect  of  receiver  field  of  view  is 
included.  For  a receiver  field  of  view  whose  full  cone  angle,  9 , exceeds  the 
full  cone  angle  of  the  clear-weather  blur  circle,  all  of  the  unscattered  light 
intercepted  by  the  receiver  pupil  reaches  the  photodetector.  On  the  other  hand, 
if  9r  ^ tr/4,  only  a fraction  sin2(9r/2)/sin2(7r/8)  of  the  diffusion-regime  light 
collected  by  the  receiver  pupil  reaches  the  photodetector.  In  Fig.  10  we  have 
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(7) 


plotted  and 

Tr^  s Tr  sin"'(0r/2)/sinZ(7r/8)I 

versus  optical  thickness  for  various  values  of  (cos  0)  using  A , 0 , L as 

r t 

before  and  0^  = 20  mRad.  We  see  from  Fig.  10  that  the  crossover  optical 
thickness,  at  which  Bucher's  scaling  laws  predict  that  the  diffusion-regime 
radiation  incident  on  the  photodetector  will  equal  the  unscattered  radiation 
incident  on  the  photodetector,  lies  between  18  and  20  for  the  assumed  values 
of  A^,  0t,  L,  0 , and  0.  0 < (cos  0)  < 0.  8.  Once  again  for  sufficiently  large 
( cos  0>  values  (e.  g.  , ( cos  0)  = 0.  875  in  Fig.  10)  there  will  not  be  a cross- 
over optical  thickness.  For  such  (cos  0)  values  there  is  little  we  can  say 
about  transmission-fraction  behavior  in  the  gap  between  the  asymptotes. 

We  shall,  for  discussion  purposes  in  the  sequel,  assume  that  a crossover 
optical  thickness  does  exist  for  our  experiment,  i.  e.  , (cos  0)  ^ 0.83. 

This  assumption  need  not,  however,  be  physically  valid.  In  particular,  using 
Deirmendjian's  scattering  functions  [26]  for  0.  7 p wavelength  propagation 
through  water  hazes  H and  M and  water  cloud  C.  1,  we  have  computed  (cos  0) 
values  of  0.  74,  0.79,  and  0.  84,  respectively.  Furthermore,  Bucher's  calcu- 
lation [25]  of  Mie  scattering  by  water  spheres  (refractive  index  1.  34)  shows 
that  0.8  <(cos  0)  $ 0.  875  for  Mie  parameters  between  30  and  200  such  as  are 
known  to  occur  (for  0.7  p wavelength  radiation)  in  fog  [27]. 

To  conclude  the  a priori  link  analysis,  we  turn  our  attention  to  the  multi- 
path  spread  predicted  for  the  supposed  crossover  optical  thickness.  The 
RMS  time  dispersion  of  the  diffusion-regime  radiation  reaching  the  receiver 
entrance-pupil  plane,  ay  computed  from  Eq.  (3)  for  r = 18,  L = 13.6  km  and 
various  (cos  0)  values  appears  in  Table  1. 

Table  1.  Diffusion-Regime  Multipath  Predictions 


( cos  0) 

0.0 

0.  6 

0.  7 

0.  8 

(psec) 

310 

146 

115 

83 

This  calculation  is  based  on  Eq.  (3)  with  t = 18,  L = 
13.6  km. 
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C.  Propagation- Link  Measurements 


The  preceding  transmission-fraction  analysis  suggests  that  the  following 
behavior  should  be  encountered  in  our  ruby-link  measurements.  Under  clear 
to  modest-visibility  conditions,  the  detected  boresight  ruby  signal  should  be 
dominated  by  unscattered  radiation.  We  should  thus  observe  an  exponential 
extinction  with  decreasing  visibility,  a sharply  peaked  angular  spectrum 
(perhaps  0.  1 mRad  full  cone  angle),  and  negligible  lime  dispersion.  As  the 
visibility  decreases,  the  detected  boresight  ruby  signal  should  eventually  be 
dominated  by  the  diffusion-regime  light.  Assuming  our  receiver  is  sufficiently 
sensitive  to  make  measurements  in  this  regime,  we  should  observe  an  alge- 
braic decay  of  transmission  with  further  decreases  in  visibility,  an  angular 
spectrum  whose  full  cone  angle  is  order  of  magnitude  1 Rad,  and  time  dis- 
persion (multipath  spread)  well  in  excess  of  1 psec. 

We  do  not  know  from  the  transmission-fraction  analysis  whether  the 
transition  from  the  unscattered  to  the  diffusion  behavior  will  occur  abruptly 
or  gradually.  It  may  be  that  the  boresight  detected  signal  follows  the  unscat- 
tered transmission-fraction  curve  all  the  way  to  the  crossover  optical  thick- 
ness, at  which  point  an  immediate  changeover  to  the  diffusion  regime  occurs. 
On  the  other  hand,  it  is  possible  that,  due  at  first  to  single  scatter  and  later 
to  low-order  multiple  scatter,  the  measured  boresight  transmission  fraction 
gradually  changes  from  the  unscattered  result  (at  very  low  optical  thicknesses) 
to  the  diffusion-regime  result  (at  very  high  optical  thicknesses).  In  the  former 
case,  we  should  expect  to  see  little  or  no  angular  spread  or  multipath  spread 
as  the  visibility  decreases  until  the  crossover  thickness  is  reached,  at  which 
point  very  large  angular  spread  and  multipath  spread  are  suddenly  apparent. 

In  the  latter  case,  we  should  expect  to  see  gradual  increases  in  both  the 
angular  and  multipath  spread  of  the  boresight  signal  with  decreasing  visibility, 
i.  e.  , no  dramatic  changes  in  the  received  signal  will  be  manifested  as  the 
crossover  optical  thickness  is  exceeded. 

The  data  we  have  obtained  are  somewhat  equivocal  with  regard  to  the 
above  possibilities.  The  observed  multipath  spread  was  never  more  than 
50-350  nsec,  and  angular  spectra  were  at  most  a few  mRad  full  cone  angle, 
in  concert  with  the  hypothesis  that  our  measurements  were  dominated  by 
unscattered  radiation.  However,  there  was  a tendency  for  the  width  of  the 
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observed  angular  spectrum  to  increase  with  decreasing  visibility.  Further- 
more, the  near-absence  of  multipath  and  the  narrow  angular  spectrum  per- 
sisted when  measurements  were  made  with  the  receiver  purposefully  misaligned 
to  reduce  the  level  of  received  unscattered  radiation.  These  observations 
suggest  a gradual  transition  from  unscattered  transmission  to  the  diffusion 
regime,  such  as  might  be  due  to  multiple  forward  scatter  (cf.  [25]).  Finally, 
it  appears  that  our  links  are  not  sufficiently  sensitive  to  reach  the  crossover 
optical  thickness  calculated  in  section  B.  2.  Thus  it  is  not  too  surprising  that 
very  large  angular  and  multipath  spreads  were  not  observed,  i.  e.  , when  the 
visibility  was  low  enough  to  put  our  links  into  the  diffusion  regime,  we  were 
unable  to  measure  any  signal. 

The  details  of  the  experiments  are  presented  below;  in  section  C.  1 we 
describe  the  data,  in  section  C.  2 we  compare  the  experimental  results  with 
the  transmission-fraction  analysis  of  section  B.  2. 

1 . Multipath  and  Angular  Spread 

Under  clear-weather  conditions  (visibility  10  miles  or  greater),  the  ruby 
laser  pulse  appears  at  the  receiver  to  be  of  150  nsec  duration  (see  Fig.  11). 

We  believe  that  the  transmitter  produces  a 70  nsec  pulse;  a pulse  monitor  in 
Cambridge  shows  a 100  nsec  duration  limited  by  the  monitor  electronics.  We 
ascribe  the  excess  pulse  duration  observed  in  clear  weather  at  Waltham  to  the 
finite  bandwidth  of  the  post-detection  electronics;  a single-photon  event  has 
a rise-fall  time  of  30  nsec.  In  low-visibility  weather  (visibility  <6  miles), 
three  types  of  received-pulse  behavior  have  been  observed  on  the  ruby  link: 

i)  pulse  shape  and  duration  identical  with  that  seen  in  clear  weather; 

ii)  pulse  duration  comparable  to  that  seen  in  clear  weather,  but  with 
internal  structure  (see  Fig.  12); 

iii)  pulse  duration  50-350  nsec  longer  than  that  seen  in  clear  weather; 
pulse  has  a sharp  rising  edge  and  a slow  decay  (see  Fig.  13). 

In  Table  2 we  have  summarized  our  multipath  observations,  using  the 
designations  i),  ii),  and  iii)  to  categorize  pulse  behavior  as  described  above. 
These  observations  were  made  with  the  transmitter  and  receiver  in  boresight 
alignment.  Measurements  taken  on  or  before  7/21/76  were  made  using  a 
1. 8 mRad  FWHM  transmitter  beam  divergence;  measurements  taken  after 
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SIGNAL  VOLTAGE  _ SIGNAL  VOLTAGE 


0.0  0.1  0.3  0.5 

TIME  (fisec) 

Fig.  11.  Clear-weather  ruby  pulse. 


Fig.  12.  Type  ii  ruby  pulse. 


Fig.  13.  Type  iii  ruby  pulse. 
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Date  I Visibility  Pulse  Characteristics  Weather  Conditions 
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11/10/76  7 miles  5-6  miles  ii  150  nsec 


7/2l/7t>  were  made  using  a 10  mltad  FWII1VI  beam  divergence.  On  all  days, 
measurements  were  made  with  various  receiver  fields  of  view  ranging  from 
1 mRad  (or  less)  to  20  mRad  full  cone  angle;  on  no  occasion  was  a variation 
in  time  dispersion  with  field  of  view  noted.  On  several  days,  measurements 
were  made  with  the  receiver  aimed  above  the  transmitter;  no  significant 
increase  in  time  dispersion  was  found,  not-  was  there  any  discernible  dif- 
ferential time  delay  between  on-axis  and  ofl-axis  observations.  (Because  of 
synchronization  .jitter,  differential  delays  of  less  than  100  nsec  cannot  be 
distinguished  from  timing  errors.) 

The  near-total  absence  of  multipath  in  the  on -axis  data  is  consistent  with 
the  hypothesis  that  these  measurements  are  primarily  sensing  unscattered 
radiation.  However,  the  absence  of  multipath  or  differential  delay  in  the  off- 
axis  data  indicates  that  the  scattered  radiation  is  itself  largely  free  of  time 
dispersion.  It  is  of  interest,  therefore,  to  examine  the  angular  spectrum 
measurements.  An  on-axis  signal  dominated  by  unscattered  radiation  will 
manifest  a sharply  peaked  angular  spectrum.  An  off-axis  scattered  signal 
that  is  essentially  free  of  multipath  and  differential  delay  should  also  exhibit 
a narrow  angular  spectrum,  for  it  would  otherwise  be  likely  for  there  to  be 
excess  path  lengths  (relative  to  on-axis  unscattered  transmission)  of  more 
than  the  100  feet  needed  to  be  resolvable  as  100  nsec  dispersion  or  delay.  We 
shall  examine  data  from  both  the  ruby  and  holmium  links. 

Under  clear-weather  conditions  (visibility  10  miles  or  greater),  both  the 
ruby  and  holmium  wavelength  angular  spectra  are  quite  narrow;  typical 
examples  appear  in  Figs.  14  and  15.  We  believe  that  in  clear  weather  the 
width  of  the  ruby  wavelength  angular  spectrum  is  due  to  atmospheric  turbu- 
lence [12],  (Indeed,  in  a stable  8 mile  visibility  haze,  narrower  ruby  angular 
spectra  are  ebserved  than  in,  presumably  more  turbulent,  12  mile  visibility 
clear  windy  weather.)  We  believe  that  the  width  of  the  clear  weather  holmium 
angular  spectrum  is  due  to  telescope  aberrations.  (The  Zoomar  optics,  and 
in  particular,  the  Schmidt  corrector  plate,  were  not  designed  for  use  in  the 
infrared.) 

Under  low-visibility  conditions,  the  ruby  and  holmium  angular  spectra 
continue  to  have  the  general  characteristics  of  Figs.  14,  15;  typical  examples 
appear  in  Figs.  16-19.  These  low-visibility  observations  are  for  coaxial 
transmitter-receiver  alignment;  thus  they  may  be  indicating  nothing  more 
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12  mile  visibility,  10/l8/7f 
IS  sec  of  boxcar  averaging 
with  4 sec  effective  time 
constant 
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2-4  mile  visibility  in  fog, 
11/5/76,  IS  see  of  boxcar 
averaging  with  4 see  effee- 
tivr  time  constant 
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than  unscattered  radiation.  To  a certain  extent  this  conclusion  is  negated  by 
the  tendency  of  the  angular  spectra  to  broaden  somewhat  in  the  lowest  visibil- 
ities Moreover,  angular  spectrum  measurements -made  witirthr  transmitter- 
purposefully  misaligned  (to  reduce  unscattered  transmission  to  the  photodetec- 
tor) show  the  same  general  character  as  Figs.  16-19,  cf.  Fig.  20.  Thus  it 
would  appear  the  scattered  radiation  intercepted  by  the  receiver  arrives  with 
a narrow  (few  mRad  cone  angle)  angular  spectrum. 


2.  Transmission-Fraction  Comparison 

We  shall  conclude  our  propagation-measurement  discussion  with  an 
attempt  to  reconcile  the  data  of  section  C.  1 with  the  transmission-fraction 
analysis  of  section  B.  2.  In  Tables  3-5  we  have  assembled  our  transmission 
data  at  ruby  and  holmium  wavelengths.  Using  the  relationships  [4] 

t = 26.  7/V  at  694.  3 nm,  (8) 

and 

t = 1 6.  8/v  at  2.  06  p,  (9) 

to  convert  visibility,  V,  measured  in  miles  to  optical  thickness  at  ruby  and 
holmium  wavelengths,  we  obtain  the  data  points  shown  in  the  transmission 
versus  t profiles.  Figs.  21-23. 

The  rubv-wavelength  data  (Table  3,  Fig.  21)  were  taken  with  coaxial 
transmitter-receiver  alignment,  10  mRad  transmitter  beam  divergence,  and 
20  mRad  cone  angle  receiver  field  of  view.  Thus,  the  experimental  ruby- 
wavelength  transmission  profile  can  be  directly  compared  with  the  unscattered 
and  diffusion-regime  profiles,  T^  and  Trf,  from  Fig.  10.  A cursory  exami- 
nation of  Figs.  10  and  21  reveals  two  immediate  features  of  the  experimental 
profile:  no  measurements  were  obtained  at  optical  thicknesses  well  in  excess 

of  the  predicted  crossover  point;  with  the  exception  of  the  two  lowest  visibility 
observations,  the  experimental  profile  seems  to  follow  a Bo uguer- Lambert 
law.  The  reason  we  have  no  data  well  beyond  the  predicted  crossover  thickness 
is  unequivocally  lack  of  system  sensitivity;  we  have  never  been  able  to  see 
the  ruby  signal  in  1 mile  or  less  visibility  at  both  terminals.  To  explore  the 
contention  that  the  experimental  transmission  profile  is  explainable  as  an 
exponential  extinction  of  unscattered  light,  we  have  least-squares  fitted  two 
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Table  3.  Ruby  Laser  Transmission  Data 
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Table  4.  Holmium  Laser  Transmission  Data:  1 mRad  FOV 
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Table  5.  Holmium  Laser  Transmission  Data:  20  mRad  FOV 
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SIGNAL  VOLTAGE 


OPTICAL  THICKNESS,  t 


Fig.  22.  Holmium-link  signal  vs  optical  thickness  (based  on  Table  4). 
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Fig.  23.  Holmium-link  signal  vs  optical  thickness  (based  on  Table  5). 
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exponential  laws  to  the  data.  The  dotted  line  in  Fig.  21, 


•y 

v = 7.  09  X 10"  exp(-0.  4 6 t), 


(10) 


is  a least-squares  fit  of  an  exponential  curve  to  all  11  data  points;  the  coef- 
ficient of  determinism  for  Eq.  (10)  is 
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where  {(vn>Tn):  l^n^  ll}  is  the  data.  The  solid  line  in  Fig.  21, 


v = ‘ . 08  X 10  exp(- 1 . 01  t), 


(12) 


is  a least-squares  fit  of  an  exponential  curve  to  the  9 data  points  with  the 

2 

lowest  t values;  its  coefficient  of  determinism,  r , equals  0.  827. 

On  the  basis  of  Eq.  (12),  we  believe  our  20  mRad  field-of-view  ruby- 
wavelength  transmission  profile  definitely  follows  the  exp(--r)  extinction  law 
for  unscatte.'ed  radiation  at  least  to  optical  thickness  10.  The  two  data  points 
we  have  at  higher  optical  thicknesses  lie  well  above  the  exp(-T)  line  at  their 
nominal  optical  thickness  values.  This  data  may  be  indicating  a transition  to 
diffusion-regime  transmission  behavior.  However,  the  enormous  error  bars 
on  optical  thickness  for  these  two  points  (which  are  due  to  the  highly  inhomo- 
geneous character  of  the  propagation  conditions  that  were  extant  11/5/76) 
permit  these  points  to  be  placed  in  close  agreement  with  Eq.  (12).  Thus,  we 
cannot  say  with  certainty  that  these  points  are  not  explainable  by  an  exponential 
extinction  lav/. 

In  Figs.  22  and  2 3 we  have  plotted  holmium-wavelength  transmission  data 
for  coaxial  tiansmitter-receiver  alignment  with  1 mRad  (Fig.  22)  and  20  mRad 
(Fig.  2 3)  cone-angle  receiver  fields  of  view.  Both  of  these  figures  appear  to 
follow  exponential  extinction  laws.  Moreover,  they  indicate,  when  compared 
with  Fig.  21,  the  decreased  optical  thickness  measured  at  2.  06  p relative  to 
that  measurei  at  694.  3 nm  (cf.  ruby  points  1,  2,  4,  9 and  holmium  points 
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4,  6,  7,  9).  A least-squares  fit  of  an  exponential  law  to  all  9 holmium  data 
points  yields  the  relationships 

v = 1.  11  X 101  exp(-0.  66  t),  (13) 

for  the  1 mRad  field  of  view,  and 

v = 1.  24  X 101  exp(-0.  63  t),  (14) 

for  the  20  mRad  field  of  view.  These  equations  (which  are  shown  as  dotted 
lines  in  Figs.  22  and  23)  have  coefficients  of  determinism  0.  728  and  0.  7 30, 
respectively.  If  we  omit  the  data  obtained  on  10/6/76,  because  of  its  sub- 
stantial error  bars  on  optical  thickness,  we  find  the  least-squares  exponential 
curves  (the  solid  lines  in  Figs.  22  and  2 3) 

v = 1.  87  X 101  exp(-0.  90  t),  (15)' 

for  the  1 mRad  field  of  view,  and 

v = 2.  05  X 10*  exp(-Q.  86  t)  (16) 

for  the  20  mRad  field  of  view,  with  coefficients  of  determinism  0.  80  3 and 
0.  809,  respectively.  It  might  seem  that  there  is  a slightly  slower  than  Bouguer- 
Lambert  decay  with  optical  thickness  to  our  holmium  data.  Because  of  the 
small  optical  thicknesses  in  our  data,  however,  it  is  more  likely  that  Eq.  (9), 
for  obtaining  2.  06  ^ optical  thickness  from  visibility,  is  in  error. 
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CHAPTER  III.  A THEORETICAL  APPROACH 


The  angular  spectra  pr  esented  in  Chapter  II  are  at  most  a few  milliradians 
wide,  even  for  measurements  made  at  about  ten  extinction  lengths  (t  = 10). 
rhis  narrowness  could  result  from  the  field  being  dominated  by  unscattered 
light,  but  that  explanation  does  not  appear  to  be  fully  consistent  with  the  data 
obtained  when  either  the  transmitter  or  the  receiver  are  pointed  off  the  optical 
axis  (see  Chapter  II). 

lire  observed  narrowness  and  the  uncertainty  as  to  whether  it  was  asso- 
ciated with  the  unscattered  field  led  us  to  develop  a theoretical  approach  for 
analyzing  multiple-scatter  propagation  in  the  region  where  the  fields  possess 
such  spectra.  That  approach,  developed  from  the  transport  equation,  is  pre- 
sented in  this  chapter.  It  differs  from  the  approach  that  has  been  used  so  suc- 
cessfully to  describe  propagation  in  turbulence  in  that  the  particle  scattering 
function  is  assumed  to  be  either  wide  or  smooth  compared  to  the  angular  spec- 
trum of  the  field. 

The  order  of  the  chapter  is  as  follows.  In  section  A the  description  of 
photon  propagation  by  the  linear  transport  equation  is  briefly  reviewed  and  the 
interpretation  of  the  solution  as  a channel  transfer  function  is  given.  In  sec- 
tion B the  equation  is  simplified  by  assuming  that  the  angular  spectrum  of  the 
field  is  narrow  relative  to  the  particle  scattering  function  and  the  solution  is 
reduced  to  the  evaluation  of  inverse  Fourier  Transforms.  Finally,  in  sec- 
tion C,  the  solution  is  obtained  in  detail  for  the  isotropic  scattering  problem. 
This  problem  is  probably  not  the  one  that  best  describes  the  scattering  environ 
ment  encountered  in  our  measurements,  but  it  does  illustrate,  in  a relatively 
simple  context,  the  approach  that  can  be  applied  to  that  environment. 

A.  The  Transport  Equation 

Our  approach  to  multiple -scatter  propagation  is  through  the  linear  trans- 
port equation  [28-30 J.  That  is,  the  "field"  is  viewed  as  a collection  of  photons 
oat  n of  which  is  characterized  by  its  position  r and  direction  of  propagation  ft 
as  functions  of  time  t. 
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1.  Formulation 


The  trajectory  in  r,12  space  followed  by  each  photon  as  time  evolves  is 
realistically  assumed  to  be  a Markov  process.  That  is,  the  statistical  behavior 
of  the  trajectory  for  t greater  than  t ' , given  knowledge  of  the  entire  trajectory 
prior  to  t ' , is  identical  to  the  behavior  given  only  knowledge  of  the  values  of 
r and  12  at  t'.  Thus  the  probability  density,  p(r,12,  t),  for  a photon  to  be  at  the 
position  r with  a direction  vector  $2  at  the  time  t satisfies  the  Chapman 
Kolmogorov  equation  [31,  32 J.  Since  changes  in  take  place  only  at  the  dis- 
crete locations  where  a photon  encounters  a scattering  particle,  this  equation 
reduces  to  the  linear  transport  equation  [28-30 J 

(c  +ilr+“'v)  p(^,“,  t)  = f ps(n-n')  pfr.fl'.t)  dn\  (17) 

In  Eq.  (17)  c is  the  velocity  of  light,  V is  the  gradient  operator  with 
respect  to  r,  S(12-12')  is  the  particle  scattering  function  and  a and  (3  are, 
respectively,  the  local  extinction  and  scattering  coefficients.  In  general  both 
c and  p can  be  functions  of  position  but  for  simplicity  we  take  them  to  be  con- 
stant throughout  the  medium.  Thus  exp  -a r is  the  probability  that  the  photon 
travels  at  least  a distance  r without  encountering  a particle,  p/a  is  the  prob- 
ability that  the  photon  will  be  scattered  upon  encountering  a particle  and 
(a-p)/c  is  the  probability  that  it  will  be  absorbed. 

The  directional  characteristics  of  the  scattering  particles  are  described  by 
the  scattering  function  S(12’12')  which  is  the  conditional  probability  density  for 
a photon  to  be  scattered  in  the  direction  12  given  it  is  incident  upon  the  scat- 
terer  with  direction  vector  12'.  Thus 

J S(n-&)  d 12  = 1,  (18) 

where  d!2  denotes  a differential  surface  element  of  the  unit  sphere.  In  spher- 
ical coordinates  with  12'  as  the  zenith  direction 

f02n  dp  dG  sin  6 S(cos  0)  = 1 (19) 

In  general,  the  scattering  function  can  depend  upon  both  12  and  12',  rather 
than  just  upon  their  inner  product  and  it  can  also  depend  upon  position.  The 
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assumption  that  it  depends  only  upon  • i2'  is  often  realistic;  the  assumption  that 
it  is  independent  of  position  is  made  for  simplicity. 

For  our  purposes  it  is  convenient  to  express  Eq.  (17)  in  terms  of  the 
I.aplace  transform,  from  t to  s,  of  p(r,Q,  t)  as  an  integral  equation  rather  than 
an  integro-differential  equation.  For  photons  produced  at  time  t = 0 with 
an  initial  probability  density  in  position  and  an  angle  of  P (r, fl)  the  result  is 
easily  shown  to  be. 

P(r, £2,  s)  = p dv  / d£2«  S(£2’£2')  P(r-v£2,£2«,  s)  exp(-vy)  + S (r,S2)  (20a) 

with 

Su(r,tt)  s fQX  dv  Po(r-vtl,i7)  exp(-yv)  (20b) 

and 

Y = a + (s/c).  (20c) 

Here  P(r,£2,  s)  is  the  Laplace  transform  of  p(r,  £2,  t); 

P(r,  s)  = Jj°  dt  p(r,£2,  t)  exp(-st).  (21) 

2.  The  Channel  Frequency  Response 

In  principle,  the  solution  to  Eq.  (20a)  can  be  inverse-transformed  to  obtain 
p(r,  £2,  t)  from  which  the  multipath  behavior  of  the  channel  can  be  inferred.  How- 
ever, it  is  simpler,  and  almost  equally  enlightening,  to  interpret  the  behavior 
in  the  s domain. 

To  develop  that  interpretation,  suppose  that  the  source  produces  photons  at 

an  instantaneous  average  rate  of  \(t)  photons  per  second  and  that  each  generated 

photon  is  distributed  in  space  and  angle  according  to  the  probability  density 

P (r,  ft).  Note  that 
o 

/ dr  / d£2  P (r.fi)  = 1 (22) 

and  that  \(t)  is  nonnegative.  We  seek  \(r,  f2,  t),  the  instantaneous  average  den- 
sity of  photons  per  second  appearing  at  the  position  r with  a direction  £2  at 
the  time  t. 

If  X.(t)  were  an  impulse  6(t--r)  X.(t)  of  area  X.(t)  at  time  t,  X.(r,£2,  t)  would 
be  \(t)  p(r, £2,  t' -t)  where  p(r, ft,  * ) is  the  inverse  transform  of  P(r, C2,  s). 
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For  the  general  time  input  \(t), 

Mr,  ft,  t)  = / dT  \(t)  p(r,  ft,  t - t)  (23a) 

or,  in  the  transform  domain, 

A(r,  ft,  s)  = A(s)  P(r,  ft,  s)  (23b) 

where  A(s)  and  A(r,ft,  s)  are,  respectively,  the  Laplace  transforms  of  X.(t) 
and  X(r,f2,t).  Thus  P(r,ft,  s),  the  solution  of  Eq.  (20),  can  be  inter- 
preted as  the  transfer  function  in  the  s domain,  between  the  time  varia- 
tion in  the  average  generation  rate  of  the  source  and  the  average  density  of 
photons  per  second  at  r with  direction  ft. 

Note  that  the  interpretation  of  P(r,  ft,  s)  as  the  transfer  function  for  the  aver- 
age photon  rate  does  not  require  that  different  "photons"  are  independent.  To 
the  extent  that  they  are  independent,  P(r,  ft,  s)  also  determines  the  higher  order 
moments  of  the  flux.  In  particular,  the  variance  of  the  local  density  of  photons 
at  r,  ft,  t will  be 

<r2( r,  ft,  t)  = Mr,  ft,  t).  (24) 

More  generally,  in  the  low  photon  coherence  approximation,  the  number  of 
photons  in  the  vicinity  of  r,  ft  at  the  time  t will  be  Poisson-distributed  with 
a rate  parameter  Mr,  t)  which  is  determined  by  \(t)  through  the  transfer  func- 
tion P(r,  ft,  s; . 

It  should  also  be  noted,  for  future  reference,  that  P(r,ft,  0)  is  just  the  time 
integral  of  pvr,ft,  t).  Thus,  from  Eq.  (23a),  it  is  the  ratio  of  the  total  average 
number  of  photons  that  are  ever  at  the  point  r with  a direction  vector  ft  to 
the  total  average  number  of  photons  emitted  by  the  source.  If  these  numbers 
are  infinite,  but  their  time  averages  are  finite,  P(r,ft,  0)  is  the  transfer  func- 
tion between  these  time-averaged  quantities. 

The  transport  equation  is  deceptively  simple  in  appearance  but  has  been 
solved  for  only  isolated  situations.  Thus,  to  gain  some  insight  into  the  experi- 
mental results  discussed  in  section  B,  the  approximation  introduced  at  the 
beginning  of  this  chapter  will  be  developed. 
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H.  An  Approximation  to  the  Transport  Equation 


upproxt null  ion  to  In  employed  was  suggest  id  by  Ihe  obsei  veil  nar- 
i owiU'.sa  o I (lit*  v\  1 1 1 1 1 1 o i p(  r,  £2,  l ) in  '•  - uml  will  In-  i it  * vi  • lopi  *i  I in  1 1 ui  | w'u  v A 
ill  1 1 1*1  1*1  II  1 on.  Ill  mu  will.  1 1 I . ■ : I . I ; . Ill  III.*  Mil  1 1 II*  i 1*1111  II  Will  I III  II  I..  iI.  mci  lb.  d. 

l.  A Narrow  Spectrum  Approximation 

Suppose  that,  in  Eq.  (20),  P(r-v£2>£21,  s)  is  concentrated  in  £2'  about  some 

direction  £2_  which  may  depend  upon  r - v£2,  denoted  by  fT.  Specifically, 

P 

suppose  that  the  integrand  in  that  equation  vanishes  when  £2'  and  £2  differ  by 

p 

an  angle  that  is  small  with  respect  to  the  angular  width  of  the  particle  scat- 
tering function  S(  • ),  To  the  extent  that  this  is  true,  P(r-  v£2,£2',  s)  behaves 
as  an  impulse  in  £2'  relative  to  S(£2  * £2* ) and  the  equation  may  be  written  as 

P(r, £2,  s)  = p dv  S(£2-£2J  P(?-v£2,  s)  exp(-vy)  + S (r,£2)  (25a) 

P U 

with 

P(r,  s)  = / P(r, £2,  s)  d£2  (25b) 

and,  as  before, 

Su(r,T2)  = /0°°  dv  PQ(r- v£2,£2)  exp(-yv)  (25c) 

Note  that  P(r,  s)  is  the  Laplace  transform  of  p(r,  t),  the  probability  density 
that  a photon  generated  at  t = 0 will  be  at  the  point  ? at  the  time  t without 
regard  to  its  direction  of  propagation.  Equivalently  P(r,  s)  is  the  transfer 
function  relating  X.(t)  to  the  average  number  of  photons  per  unit  volume  at  the 
position  r at  the  time  t without  regard  to  their  direction  of  propagation. 

It  is  instructive  to  note  that  P(r,  £2,  s)  will  be  narrow  with  respect  to  £2, 
as  claimed,  only  if  the  integral  in  Eq.  (25a)  vanishes  whenever  £2  differs  much 
from  £2^;,  the  direction  about  which  P(r,  £2,  s)  is  assumed  to  be  concentrated. 

In  particular,  P(r,  £2,  0),  the  transfer  function  for  average  power  or  energy, 
will  be  narrow  only  if  the  nonnegative  integrand  satisfies 

S( £2  - £2__)  P(r  - v£2,  0)  exp(-va)  « 0 all  v (26) 

P 
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for  any  if  that  differs  from  if  by  more  than  the  nominal  half-width  0 of 
— ^ r c* 

P(r,  if,  0). 

For  large  enough  values  of  v the  required  vanishing  results  from  the 
exponential  factor  va.  For  smaller  values  of  v,  m particular,  for  values  of 

av  less  than,  say,  five,  either  S(£f-ifJ  or  P(?-vtT,  0)  must  vanish  whenever 

- p 

if  differs  from  if _ by  more  than  0 . Since  this  width  is  assumed  to  be  sub- 
r c 

stantially  less  than  that  of  S(  • ),  the  required  vanishing  can  be  caused  by 

S(  if  -if_J  only  if  the  angle  between  if _ and  if  nearly  exceeds  the  halfwidth  0 
P ^ P S 

of  S(  • ).  Otherwise  it  must  be  P(r-vif,  0)  that  vanishes.  These  constraints  are 

summarized  in  Fig.  24. 

The  requirements  of  Fig.  24  can  be  satisfied  in  a variety  of  ways,  depending 

upon  the  geometry  of  the  source  and  observation  region.  For  the  problems 

of  interest  here  r is  typically  located  near  the  axis  of  symmetry  of  the  source 

angular  pattern.  Were  it  not  for  scattering  the  position  r would  be  illuminated 

by  the  transmitted  beam.  For  this  situation,  Fig.  24  requires  that  any  "beam" 

that  continues  to  exist  in  the  presence  of  scattering  must  be  such  that  the 

portions  of  it  which  lie  outside  the  cone  of  half-angle  Gc  and  within  a distance 

of  a few  a of  r must  be  confined  to  the  region  where  ff_  and  if  differ  by  more 

than  0 . P 

s 

For  this  requirement  to  be  met  exactly  it  is  necessary  that  0 exceed 
0g-  in  violation  of  the  assumption  that  led  to  Eq.  (25).  However,  it  will  be 
approximately  satisfied  if  the  beam  diameter  is  small  with  respect  to  the  cone 
shown  in  Fig.  24,  i.  e.  , if  the  beam  radius  at  r is  less  than  a few  times 

a tan  0 . 

c 

The  foregoing  argument  suggests  that  the  angular  width  of  P(r,if,  s)  will 
be  narrow  as  assumed  only  if  the  beam  diameter  is  relatively  small.  This 
might  in  turn  suggest  that  Eq.  (25),  which  results  from  that  approximation, 
is  essentially  a single  scatter  approximation.  In  fact,  it  is  not.  Indeed 
Eq.  (25)  is  exactly  the  same  as  Eq.  (20)  for  isotropic  scatter;  that  is,  when 
S(if 'if')  equals  l/(47T)  for  all  values  of  if.  Thus  there  is  some  reason  to  hope 
that  Eq.  (25)  may  provide  a useful  approximation  to  P(r,iT,  s)  even  when  its 
angular  width  is  not  narrow. 

That  hope  is  reinforced  by  the  fact  that  Eq.  (25)  can  also  be  obtained  from 
Eq.  (20)  when  S(  • ) is  sufficiently  smooth.  Specifically,  the  integral  in  Eq.  (20) 


49 


can  be  evaluated  by  parts  to  yield  an  equation  which  differs  from  Eq.  (25)  only 
through  the  presence  of  an  additional  integral  whose  integrand  involves  products 
of  the  partial  integral  of  P(r,it,  s)  with  respect  to  SI  and  the  partial  derivatives 
of  S(n*S2')  with  respect  to  the  components  of  Q' . To  the  extent  that  the  deriva- 
tives are  negligible,  Eq.  (25)  is  obtained. 


2,  The  Solution  for  P(r,  s)  and  P(r,S2,  s) 

The  principal  advantage  of  Eq.  (25)  relative  to  Eq.  (20)  is  that  it  can  be 
solved  more  simply.  In  particular,  a Fredholm  integral  equation  for  P(r,  s) 
can  be  obtained  by  integrating  Eq.  (25)  with  respect  to  SI.  The  result  is,  upon 
introducing  p“  = r - vfl  and  shifting  to  a volume  integration. 


P(r,  s)  = / dp  K(r  - p* ) P(p,  s)  + Su(r) 


with 


K(r-p)  = (3  S(w) 


exp  -yj  r-p 


l-  -.2 

r-p 


(27) 


(28) 


(r-p)  _ 

w = 

(f-p|  p 

and 

r dp  ^ _ 

su(r)  = j — exp(-yp)  PQ(r-p,p/p).  (29) 

Since  is  the  direction  about  which  P(p,f2,  s)  is  concentrated  in  £1,  it 

P _ 

must  be  determined  before  Eq.  (27)  can  be  solved  for  P(r,  s).  In  principle, 

£?_  can  be  determined  from  Eq.  (25),  but  for  the  beam  geometries  that  will  be 
P 

considered  here,  it  can  reasonably  be  assumed  to  be  along  the  nominal  direc- 
tion of  the  beam.  We  take  that  direction  to  be  the  positive  z direction  so  that 


fi  *1  . (30) 

P Z 

With  that  choice  Eq.  (27)  can  be  solved  by  Fourier  transform  techniques. 

Specifically,  let  &(Z,  s),  k(Z),  and  <Su(w)  be,  respectively,  the  three- 

dimensional  Fourier  transforms  of  P(r,  s),  K(r),  and  Su(r);  e.  g. , 
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. >’(u>,  s)  / dr  P(r,  s)  i'\p(-|u>*  r) 


( Ha) 


and 

P(r,  s)  = ( 2^)  \ ^(w,  s)  cxp(  )u-  r). 

In  the  transform  domain  Eq.  (27)  can  be  written  as 
~P(u>,  s)  = k(w)  ./>(w,  s)  + ^u(w). 

Thus 

,s  (00) 

.^(u>,  S)  = — 

1 - k(w) 

or,  more  conveniently  for  our  purposes, 

- *(<*>)  _ 

.^(co,  s)  - «\(w)  + ^ S (w) 

1 - k(w)  u 


(31b) 


(32) 


It  follows  immediately  from  Eq.  (32)  that  P(r,  s)  is  given  by 
P(r,  s)  = SJr)  + p / dp  h(r-pT)  Su(p) 


(33a) 


with 


/ 1 \3  r k <w)  _ 

h(F)  = (zf)  T : ^ exP(Jw‘  r) 

1 “ k(w) 


(33b) 


The  response  function  h(r)  can  be  interpreted  as  follows.  Suppose  the 
source  function,  PQ(r , S2),  equals  6(r)  S(sT - £ ),  i.  e.  , a point  source  with  an 
angular  pattern  S(f2  • 1^)  is  located  at  the  origin.  Then  Su(r)  equals  K(r)/P. 
>u(w)  equals  k(w)/B  and,  for  this  source,  P(r,s)  equals  h(r)  as  defined  by 
Eq.  (33b).  That  is,  h(r)  is  P(r,  s)  for  a point  source  located  at  the  origin 
whose  angular  distribution  is  S(S2  • £^).  It  is  as  though  the  source  is  a single 
scatterer,  located  at  the  origin,  illuminated  by  photons  propagating  along  the 
z axis. 

For  a general  source,  the  first  term  in  Eq.  (33a)  is  the  contribution  to 
P(r,  s)  of  the  photons  that  have  not  been  scattered  and  the  second  term  is  the 
integral  of  the  product  of  pS(fT),  which  is  the  probability  that  a source  photon 
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d 


is  first  scattered  at  the  point  p and  h(r-p),  which  is  the  contribution  of  that 
"point"  source  to  P(r,  s). 

Once  P(r,  s)  has  been  determined  through  Eq.  (33),  P(F,f2,  s)  can  be  eval- 
uated from  Eq.  (25). 

Time  has  not  permitted  an  exploration  of  Eqs.  (3  3)  nor  of  the  domain  of 
validity  of  Eq.  (25).  However,  the  solution  for  the  isotropic  scatter  problem 
has  been  obtained.  For  it,  we  will  find  that  the  angular  spectrum  will  not  be 
narrow  when  a r is  much  larger  than  one.  Also,  when  nr  is  substantially  larger 
than  one,  the  bandwidth  over  which  P(r,  s)  is  constant  in  magnitude  is  limited 
to  about  0.  32  cp/farF  rad/sec  when  a = p and  about  1.  8 csTp/r  [ (o/p)  - 1 
rad/sec  when  c and  p are  nearly,  but  not  exactly,  equal.  Thus,  as  might  be 
expected  from  the  scattering  characteristics  of  hazes  and  fogs  [l,9,  26 j,  our 
experimental  data  are  consistent  with  the  theoretical  results  for  isotropic  scat- 
ter only  when  the  field  is  dominated  by  unscattered  radiation. 
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C.  Isotropic  Scatter 


The  isotropic  scatter  problem,  Cor  which  S(U  • f ) equals  1/4tt,  has  been 
considered  by  others,  beginning  apparently  vvitli  Case  et  a I . [ 3 3 J . The  results 
presented  here  differ  from  those  earlier  ones  primarily  through  the  presence 
of  the  Laplace  variable  s. 


1.  Isotropic  Point  Source 
For  isotropic  scatter  Eq.  (28)  becomes 
S(£2-T  ) = 1/477 

_ P exp(-yr) 

E(r)  = 2 (34) 

477r 

_ P V + jw 

k(w)  = -r—  In : — (35al 

2 jw  y - ,iw 

or  equivalently 

k(w)  =-tan_1y.  (35b) 

Here  u denotes  the  magnitude  of  w and  the  branch  of  the  logarithm  is  chosen 
so  that  k(u>)  = P/y  for  cj  - 0. 

Introducing  Eq.  (35)  into  Eq.  (33b)  yields 


h(r) 


tan  ^w/y) 


K _ 1 

1 , 1 (jl) 

1 tan  — 

u>  Y 


exp(jw  r). 


(36) 


For  s = 0,  the  integral  in  Eq.  (36)  reduces  to  that  evaluated  by  Case 
et  al.  [3  3].  The  evaluation  for  general  imaginary  values,  |wm>  of  s is  carried 
out  in  Appendix  I.  The  result  is,  for  nonnegative  values  of  wm> 


h(r)  = 


A(«m)  exP(Jajpr) 


2irr 


+ Mflr,  jwm) 


exp(-ar) 

~ 2~ 
47Tr 


(37) 


where 
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C 2 

to 

p 

r 

Pv 

pj1 

2 , 2 

V + w 

Y P 

0 

Here, 


i'or  U>  Si  W 
m c 


v = a + — = a + j 

t c J 


is  the  unique  positive  solution  of 


- /-f^i  = Coth 

“ V2“c/  \ c 


and  Wp  is  that  unique  solution  of 
1 + ( jw  /v) 

p'  " y 

In = q- 

2(j“p/\)  1 * (JWp/v) 


(38a) 


(38b) 


whose  imaginary  part  is  nonnegative  with  the  branch  of  the  logarithm  chosen 
so  that  In  (1)  is  zero.  The  function  e(ar,  iwm)  is  given  by 

r fr  P .!  f /irp\2]_i 

€(ar,  ja>  ) = yr(exp  or)  \ dv  ^ I 1 - — - coth  v + h- 7)  ' exp(-yrv). 


C \ p -1 

€(ar,  j«m)  = yr(expcr)  dv  j 1 - — coth 


Since  h(r)  is  the  Fourier  transform  (from  t to  um)  of  a real  quantity,  it 

possesses  conjugate  symmetry  in  Thus  the  above  results  determine  its 

value  for  negative  values  of  as  well  as  positive  values. 

The  first  term  in  Eq.  (37)  does  not  tend  to  zero  as  w tends  to  u . In 

me 

fact,  as  shown  in  the  appendix,  then  tends  to  (7rPyc)/(2wJ.  Thus  it  might 
appear  that  h(r)  changes  discontinuously  at  u = w . However,  it  can  be  shown 

XY"t  C ^ ^ 

that  f(ar,  .)wm)  also  changes  discontinuously  as  u>m  changes  from  to  u>c  and  that 
h(r)  is  continuous,  as  one  would  expect. 
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Thr  term  « (or,  ju>  )[e\p(-yr)|/| 4ffr  | in  the  r\|ircssiuii  for  li(r)  depends 
u(xm  a,|3,  r anti  w in  a complicated  way.  However,  it  contributes  little  to 
h(r)  when  or  is  appreciably  larger  than  one,  so  little  that  we  will  subsequently 
neglect  it.  To  justify  that  neglect,  note  that  h(r)  evaluated  for  = 0 gives 
the  total  average  number  of  photons  that  are  ever  in  the  vicinity  of  the  point 
r when  the  source  is  an  impulse  in  time  and  gives  the  average  number  of 
photons  per  second  there  when  an  unmodulated  source  is  employed.  The 
contribution  of  the  term  in  question  to  either  of  these  averages  is 
((or,  0)[exp(-cr)]/[4n'r1'J. 

Examination  of  Eq.  (41)  shows  that 


/-*  00 

((or,  0)  [or  exp(ar)J  \ 


or,  upon  evaluating  the  integral, 


C(ar,  0) 


Ci 


(42) 


(43) 


Thus  the  average  energy  or  power  carried  by  the  term  of  Eq.  (37)  containing 
(jar,  jw^^)  decays  with  r at  h ast  as  rapidly  as  the  unscattered  field.  Moreover, 
its  value  is  bounded  by  five  firm  that  of  the  unscattered  field  if  or  exceeds 
one  and  (3/a  exceeds  2/jr.  For  the  low  absorption  problems  of  interest  to  us 
P A will  always  exceed  2/n . Thus  the  term  in  question  can  be  neglected  almost 
as  safely  as  can  the  unscattered  field.  Henceforth  it  will  be. 

To  the  extent  that  the  second  term  of  Eq.  (37)  is  neglected, 


A(u>  ) exp(-arri) 
m r 

h(r> TttY <44> 

where 

h = “jw  /«• 

In  this  approximation,  h(r)  is  zero  for  wm  > since  A(u>m)  then  equals  zero. 
For  the  behavior  of  h(r)  is  dominated  by  exp(-arr|)  when  ar  is  sub- 

stantially greater  than  one.  Note  that  since  is  that  solution  of  Eq.  (40) 
whose  imaginary  part  is  nonnegative,  the  real  part  of  q is  nonnegative. 


56 


For  equal  to  zero  q is  real.  Its  value,  as  originally  determined  by 
Case  et  al . [33]  is  plotted  as  a function  of  p/a  in  Fig.  25.  Note  ihat  T|  van- 
ishes when  a/p  equals  one,  i.  e.  , when  there  is  no  absorption.  In  that  instance 
h(r)  does  not  decay  exponentially  with  r but  only  as  l/r.  This  corresponds 
to  the  diffusion  mode  discussed  in  Chapter  II. 

When  both  absorption  and  scattering  are  present  (a  * p * 0)  h(r)  does  decay 
exponentially  with  r but  less  rapidly  than  does  the  unscattered  field.  However, 
for  the  decay  rate  to  be  no  larger  than  one  fifth  of  the  unscattered  decay  rate 
(r|<0.  2),  p/a  must  exceed  0.985.  That  is,  only  a small  amount  of  absorption 
can  be  tolerated  if  the  decay  rate  qar  is  to  be  substantially  less  than  that  of 
the  unscattered  field. 

The  variation  of  wc>  the  frequency  above  which  M^m)  vanishes,  is  shown 
in  Fig.  26  as  a function  of  p/a.  As  noted  earlier,  when  u>m  equals  u>c, 
equals  (.TPyc)/(2w  ).  Thus  at  = u>c  the  real  part  of  q equals  (7rP)/(2a)  and 
the  imaginary  part  equals  -(JTpc)/(2w^).  Note  that  u. /(Pc)  is  on  the  order  of 
1.  4 for  P/a  greater  than  0.  5 and  the  imaginary  part  of  q is  then  about  one.  For 
P nearly  equal  to  a the  real  part  of  q is  approximately  1.  6. 

In  Figs.  27  and  28  the  real  and  imaginary  parts  of  q are  plotted  as  func- 
tions of  wm/(cP)  for  a/p  = 1,  1.04,  and  1.  11.  These  ratios  correspond  to 
ratios  of  absorption  to  scattering  cross  sections,  (a-p)/p,  of  0,  0.  04,  and 
0.  11.  Note  that  as  the  absorption  increases  with  P fixed  the  real  and  imaginary 
parts  of  q vary  less  rapidly  with  w for  u;m/(Pc)  S0.1.  Thus  the  useful  band- 
width  increases  with  adsorption  as  one  expects.  However,  this  increase  is 
obtained  at  the  expense  of  an  increased  attenuation  at  all  frequencies. 

Although  it  is  not  apparent  in  the  log-linear  plot  of  Fig.  28,  the  imaginary 
part  of  q is  a nearly  linear  function  of  u for  <"m/(cP)  less  than  about  0.  04 
and  0.11  when  (a-p)/p  is,  respectively,  0.04  or  0.11.  This  can  be  seen  in  Fig.  29 
where  part  of  Fig.  28  has  been  replotted  on  a linear  scale. 

The  same  conclusion  is  reached  by  noting  that  the  asymptotic  solution  of 
Eq.  (40)  for  a « p and  w small  yields 


*9* 

The  other  factor,  ar,  in  the  exponent  will  increase  somewhat  with  the  increase 
in  absorption  just  described,  but  for  p/a  near  one  the  effect  will  not  be  signi- 
ficant. 
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1000 


(45) 


(3  A " P 

y ~p~ f j 


m 
cp  ' 


Examination  of  this  expression  sliows  that  the  real  part  of  r|  will  be  nearly 
constant  and  the  imaginary  part  nearly  linear  in  ^^/(cP)  when  w /(cp)  is 
substantially  less  than  (a-f3)/p,  i.  e.  , less  than  0.  04  and  0.  Oil  for  the  two 
situations  just  discussed. 

Equation  (45)  can  also  be  used  to  estimate  the  nominal  3 dB  bandwidth 
of  the  exponential  factor  exp(-arr|)  when  or  is  substantially  larger  than  one. 

For  then  exp(-crr|)  will  change  by  a factor  of  two  when  r|  changes  only  slightly, 
and  the  bandwidth  can  be  estimated  from 


0.  69  5 = (or)  - x/T 


(46) 


where  denotes  the  value  of  for  which  the  real  part  of  (arp)  increases 
from  its  zero  frequency  value  by  0.  693. 

When  there  is  no  absorption  (a  =p)  the  solution  of  Eq.  (46)  is 


“b  0.  32 


On  the  other  hand,  if  there  is  some  absorption  and  a r is  enough  larger  than 
one,  the  solution  approaches 

“b_  1.8  fa  (Q  ' V 

C(3  " P V P 

For  large  optical  thicknesses,  sr,  the  increase  in  bandwidth  that  results  from 
the  presence  of  a small  amount  of  absorption  is  substantial. 


2.  Ray  Source 

The  density  P(r,  s)  is  equal  to  h(r)  for  a point  source  that  radiates  uniformly 

in  all  directions.  For  a point  source  with  an  angular  radiation  pattern 

S (S2  • ft  ),  P(r,  s)  is,  by  virtue  of  Eq.  (33a), 
o z 

P(r,  s)  - Su(r)  + p / dp  h(r-p)  Su(p), 
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(IV) 


where  S ( r),  the  mi: 
u 


altered  field  is,  from  lq.  (29), 


S 

u 


(<•') 


exp(-Vi')  / r ' t \ 

o2  "X  r / 


We  ignore  the  unscattered  field  term  Su(r)  in  the  expression  for  P(r,  s) 
and  focus  on  the  remaining  term.  That  term  has  a component  involving  the 
function  €(  • , defined  in  Eq.  (41).  As  discussed  in  conjunction  with 
Eq.  (43),  its  contribution  is  not  significant  when  cr  is  large.  Thus,  as  there, 
it  will  be  neglected  here.  The  remaining  term  is 

r dP  S (p) 

P(r,  s)  - PA(«m)  \ — ■ - _ exp(jw  | r-p|  ).  (48) 

^ 2ir  I r — p | F 

The  properties  of  P(r,  s)  have  already  been  determined,  through  h(r),  for 
an  isotropic  point  source.  To  see  the  effect  of  source  collimation,  the  other 
extreme  of  a ray  source  will  be  considered. 

Specifically,  suppose  that  So(f2  * is  an  impulse  in  solid  angle,  in  the 

direction  jf  . Then 
z 

exp(-Yz) 

Su(r)  = 2 6<P)> 

z 

where  6(  • ) is  a two-dimensional  impulse,  p is  the  component  of  r perpen- 
dicular to  and  z is  the  projection  of  r on  Consequently, 

PA(«m)  poo  exp(j“p  //p2+(z-v)2) 

P*r’ s)  ~ \ dv  , _ IZ exp(-yv),  (49) 

0 ,/  2 ~2 
Vp  + (z-v) 

where  p is  the  length  of  (T. 

The  integral  cannot  be  evaluated  in  closed  form  but  its  essential  proper- 
ties can  be  displayed  through  integration  by  parts.  To  this  end,  define 

f~ 2 2 

-yv  + jo,  V P + (z-v)  = f(v) 

and  write 
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p(r.  s)  = - -2i~  3q  dv 


exp  f(v)j 


/2  2 

p + (z-v)  h |wp(z-v) 


where  f’(v)  denotes  the  derivative  of  f(v).  Integration  by  parts  yields 

(1  “(ft) 


P(r,  s)  = 


PA(«m)  exP('|u,pr) 


27T(yr+,|u)pZ) 


(50a) 


with 


(ft  = (yr  + jw  z)  \ dv 
P Jn 


where 


y(z-v)  + jup  h/ p2  +(z-v)“  j exp- 

• 

LJ"p 

/ f>2  + (z-v)2  -r 

-yvj 

*/ p2  + (z-v)2  | jw  (z-v)  + y +(z-v)2l 

L ^ J 


(50b) 


/ 2 2 
r = //  p + z . 

It  is  not  difficult  to  show  that,  asymptotically, 


yz  + jwpr 

(ft g (51) 

(yr  + j«pz) 


as  z tends  to  plus  infinity  with  either  p or  p/z  fixed  and  nonzero.  It  can 
further  be  shown  that  when  up  equals  zero,  (ft  is  a function  of  only  a z and  p 
The  behavior  of  (ft  strongly  suggests  that  P(r,  s)  approaches 


P/(27T) 


A(w  ) 
m 


Y + (jwpz/r) 


exp(jwpr)| 


rapidly  as  yr  exceeds  unity  — at  least  when  p is  not  zero.  The  behavior  for 
p = 0 is  more  complicated  sin  e P(r,  s)  contains  some  unscattered  contribution 
along  the  z axis  (p  = 0).  However,  if  this  component  is  locally  integrated  over 
a small  area  about  the  axis,  it  is  found  that  the  result  converges  to  the  local 
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integral  of  the  above  expression.  Since  tliut  expression  is  just  h(r)  divided 
by  [(y/P)  i- (ju)pz/pr)j,  the  i'iehl  differs  from  that  produced  by  an  isotropic 
point  source  only  by  [(y/p)  + (jw  z/pr)]  once  yr  appreciably  exceeds  unity. 
Indeed,  when  a = p and  = 0,  in  which  case  w = 0,  Eq.  (50a)  becomes 

P(r,  s)  ~ h(r)  1 

V av  J 


3.  The  Angular  Spectrum 

Finally,  let  us  determine  P(r,W,  s)  from  Eq.  (25).  Since  P(r,  s)  for  an 
isotropic  point  source  and  a ray  source  differ  only  slightly,  an  isotropic 
source  will  be  considered.  Also,  the  unscattered  component  of  P(r,  ft,  s)  will 
be  ignored  and  only  the  dominant  term,  A(wm)[exp(,iu>  r)]/(2jrr),  of  P(r,  s)  will 
be  used  in  the  integral.  Subject  to  these  conditions,  Eq.  (25)  becomes 

) 

PA(U  I r»  e*p<jUp|F-vn|) 

P(r,fi,s)  = — \ dv — — exp(-yv).  (52) 

Bit"  J 0 | r- vq| 


Upon  denoting  by  z the  projection  of  r on  t!  and  by  p the  length  of  the 
component  of  r perpendicular  to  f2,  Eq.  (52)  can  be  restated  as 


P(r,Q,  s) 


pA(wm) 

Bit" 


exp(jWp  fJ p2  +(z-v)2  ) 


dv 


2 , . .2 
P + (z-v) 


exp(-yv) 


The  right-hand  member  of  this  equation  differs  from  that  of  Eq.  (49)  only 
by  a factor  of  4 ir.  Thus  the  argument  that  led  from  Eq.  (49)  to  Eq.  (50)  can 
be  repeated  here.  Specifically,  if  z tends  to  infinity  with  either  p or  p/z 
fixed,  P(r,  f2,  s)  converges  to 


P(r,  £2,  s) 


PA(wm>  exP(Jwpr> 


8 n (Yr  + jojpz) 


Yz  + Jwpr 
(Yr  + jwpz)2 


(53) 


As  with  Eq.  (50b),  P(r,  ft,  s)  includes  an  impulsive  term  when  p = 0 (S2  in  the 
direction  of  r).  However,  if  P(r,  ft,  s)  is  integrated  over  a small  solid  angle 
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about  the  direction  r,  the  singularity  is  removed  and  the  result  is  given  by 
Eq.  (5  3)  integrated  over  the  same  solid  angle. 

Finally,  it  is  convenient  to  express  Eq.  (53)  in  terms  of  the  angle  0 
between  ft  and  r.  Then  z = r cos  0,  p = r sin  0,  and  the  conditions  on 
Eq.  (5  3)  become  that  r increases  with  either  0 or  r sin  0 fixed.  The  eval- 
uation itself  becomes 


P(r,  ft,  s)  ~ 


PA(w  ) e\|>(  iw  i*) 

m p 


y COS  0 + Ju> 


1 - 


47Tr(y  cos  0)  r(y  + Joj^  cos  0) 


or  equivalently 


P(r,  ft,  s)  ~ 


1 -- 


y cos  0 + |w  ~] 

• p 1 


47T (y  + jw  cos  0)  | r(y+.|t*)  cos  0) 


P(r,  s). 


(54) 


Note  that  P(r,  ft,  s)  is  a very  weak  function  of  0 whenever  jw  differs 

P 

substantially  from  y.  Indeed,  if  = 0 and  a = p,  in  which  case  = 0, 


P(r,  ft,  s) 


P(r,  s) 
4tr 


cos  0\ 
or  / 


and  P(r,  ft,  s)  quickly  becomes  uniform  in  ft  as  or  exceeds  one. 
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The  task  at  hand  is  to  evaluate  Eq.  (in)  vvilh  .s  = lu)  for  co  '0.  To 

1 4 ' m m 

this  end  it  s convenient  to  restate  that  equation  as 


exp(-yr)  i / i \3  f K M 

h(r ) = 2 — + IT  ' lliT  / ' dt0 exp(jw-r),  (1-1) 

47Tr  1 - x(co) 


where 


= 


2jco 

or  equivalently. 


Y + Jco 

In — 

V “ Jco 


Mco) 


(I-2a) 


(I-  2b) 


and  y equals  a + (jwm/c).  Here  k(w)  has  been  written  as  k (co)  since  it  depends 
on  OS  only  through  co. 

Upon  expressing  the  integral  in  polar  cooruinates  with  r as  the  zenith 
direction  and  evaluating  the  integrals  on  4>  and  0,  one  obtains 

2 

r 00  k (co)  exp(-yr) 

h(r)  : — \ dco  [co  sin  (cor)  J + — . 

(2JT)  pr  ^0  1 - x(w)  477 r 


Since  x(w)  is  defined  and  analytic  in  the  strip  |lm  wj  < Re  y,  and  since  both 
co  sin  wr  and  k(Co)  are  even  functions  of  co,  this  may  also  be  written  as 

2 

j p 00  cok  (co)  exp(-yr) 

h(r)  = n \ dco  exp(jcor)  + ^ — . (1-3) 

(277)^  .1(3 r -<*>  l-x(co)  477 r 


The  integral  can  be  further  reduced  by  contour  integration  in  the  complex 
co  plane.  To  that  end,  note  that  k(co).  continued  beyond  the  strip  [im  co|  < a,  is 
analytic  in  the  whole  upper  half-plane  except  for  a branch  point  at  co  = y and  for 
isolated  poles  at  the  zeros  of  1 - k(co).  Moreover,  cox^M/ll  -k(co)|  vanishes  as 
the  imaginary  part  of  co  goes  to  infinity.  Thus  the  integral  around  the'  contour 
in  the  cut  plane  shown  in  Fig.  1-1  equals  the  sum  of  the  integral  Ij  along  the 
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real  w axis,  the  integrals  1,  and  1^  along  the  cut  and  the  integral  I around  the 
branch  point.  Furthermore,  the  integral  1^  can  be  shown  to  vanish.  Thus  the 
integral  around  the  contour  reduces  to  + I.,  + K where  lj  is  the  integral 
appearing  in  Kq.  (1-3)  This  contour  integral  equals  2rrj  times  the  sum  of  the 
residues  of  the  integrand.  Thus 

exp(-y  r) 

h(r) = (I0  + I,)  + 27Tj  S Residue  («  ) (1-4) 

47?r  L 5 p 

where  the  sum  is  over  the  poles,  co^,  of  the  integrand  within  the  contour. 

The  integrals  and  1^  differ  only  in  that  they  are  on  opposite  sides  of  the 
cut  and  are  in  opposite  directions.  When  they  are  added  one  obtains 


-y  r>  oo 

’z  + h'WJ,  dl' 


r u t 

1 - - — coth  v 


p -i  T 

1 coth  v ( + 

yv 


exp(-yrv) 


1 2 + T3 


exp(-vr) 


y c 00  r P 

+ — r ) dv  { 1 

47rr  J,  ' vv 


coth  1 v 


exp(-y  rv) . 


'r°r  w equc  1 to  zero  the  last  integral  is  essentially  one  encountered  by  Case 
•it  al.  f 1 ].  F ollowing  them  we  define 


€(ar,ju/m)  = v r(exp  cr) 


r 00  fr  p 

\ dv  1 

J1  IL 


coth  ^ v 


exp(-y rv) 


and  write 


h + h*- 


exp(-yr)  e(ar,  jw  ) 

g — + exp(-ar). 


The  function  (far,  0)  is  identical  to  the  function  f(or)  used  by  Case  [ 1 ). 

It  remains  to  evaluate  the  residue  term  in  Kq.  (1-4).  The  poles  are  the 
upper  half-plane  zeros  ^ of  1 - k(go),  i.  e.  , they  are  the  solutions  of 
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1 


0-7) 


1 ' (Jjuip/V) 

In 

2(jo)  A)  I - (jw  /v) 


\ 


for  whirl)  Im  w ^ is  positivr  The  number  of  such  solutions  is  given  by  the 
number  of  times  the  y/p  point  in  the  complex  plane  is  encircled  by  the  locus 
of  [2j(cj/\)-1  In  { [ 1 -.i(w/y)|/(1  "(Wy)]}  as  co  varies  around  the  contour  shown 
in  Fig.  I-l.  That  locus  has  the  general  shape  shown  in  Fig.  1-2.  Thus  there 
will  be  one  upper  plane  solution  of  Eq.  (1-7)  if  y/p  falls  in  the  shaded  area 
and  none  if  it  does  not. 

There  is  always  a solution  when  com  equals  zero  since  y/p  is  then  real 
and  positive  and  hence  is  contained  in  the  shaded  area.  More  generally,  a 
solution  will  exist  for  small  values  of  wm  but  not  for  large  values.  The  tran- 
sition occurs  when  o>m  is  such  that  the  point  y/p  lies  on  the  locus  shown  in 
Fig.  1-2.  The  value  of  co  for  which  this  occurs  will  be  denoted  by  u . It  is 

the  value  of  w for  which  the  solution  <*>  of  Eq.  (1-7)  lies  on  the  left  side  of 
m p 

the  cut  in  Fig  I-l  Thus  ooc  is  the  real  nonnegative  solution  of  Eq.  (1-7)  for 
some  v greater  than  1 when  w equals  jyv.  Specifically, 


2 v 


, 1 - v c 

n 1 + V = p 


+ J 


Pc 


with  v > 1 


or 


1 . v + 1 

r — In r + 


F 


v - 1 2v 


w 

(L_  C_ 

p J Pc  ’ 


where  the  branch  of  the  logarithm  has  been  chosen  to  be  on  the  left  side  of  the 
cut.  The  imaginary  part  of  this  equation  implies  that 


- w 

JL  - c 

2v  " pc  ' 


d-8) 


Using  this  result  to  eliminate  v from  the  real  parts  of  the  equation  yields 


P V 


7Tpc\  . (rrpc)/(2w  ) - 1 
' = 4-  In 2 


(7rpc)/(2wc)  + 1 


or,  equivalently. 
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Fig.  1-2.  Locus  of  f2j(w/\)]  1 In  {[  1 + (jw/v)  |/[  1 - (jw/\) ) } for  co  on  the  contour  of  Fig. 


- (in)  = coth  ( \.  (1-9) 

a\Z*c)  \Zuvl 

This  is  tlu'  defining  equation  tor  <*>  . 

Since  the  hyperbolic  cotangent  varies  from  infinity  to  unity,  Eq.  (1-9)  will 

possess  a single  positive  solution.  Thus  there  is  a single  value,  u>c,  of  u>m  for 

which  the  y/p  point  in  Fig.  1-2  lies  on  the  boundary  of  the  shaded  region.  For 

all  smaller  values  it  lies  within  the  shaded  region,  and  1 - k(«)  possesses  a 

single  upper  half-plane  zero.  For  all  larger  values  the  y/p  point  is  outside 

the  shaded  region  and  there  is  no  upper  half-plane  zero.  We  note  in  passing 

that  when  <*>  = w the  zero  w is  located  at  (7rpyc)/(2to  ). 

m c p c 

It  remains  to  evaluate  the  residue  of  wp  that  appears  in  Eq.  (1-4).  Since 
cop  is  a simple  pole,  its  residue  is  given  by 


Up*  (u  ) exp(jWpr) 


lies  (Wp) 


(Ztt)  j(3r 


-K-(Wp) 


where  k'(w)  is  the  derivative  of  k(u)  with  respect  to  w.  Since  <c(wp)  = 1 and 


1 YP  1 

k'(oj)  « (co)  +~ s 5" » 

w w y2  + co2 


this  can  also  be  written  as 


Res  (gj  ) 
P 


Wp  exp(jwpr) 


(2tt)  jp  r ] 1 £ 2 

V +Wp 


(I-  1 0) 


Introducing  Eqs.  (T - f>)  and  (1-10)  into  Eq.  (1-4)  yields 


h(r)  = 


A(w  ) exp(-ar) 

exp(ju  r)  t-  €(ctr,ja>  ) — — ? — , 


where 
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r 


CO 


AKn>  - < 


Pv  1 
**  1 " ~2~  2 ' 
, +Wp 


for  0 « w < w 
m c 


for  o>  > co 
m c 


and  efar(jconi),  oop,  and  W(  are  given  by  Eqs.  (1-5),  (1-7),  and  (1-9). 
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METRIC  SYSTEM 


BASE  lINTrS: 


Quantity 

Unit 

SI  Symbol 

Eormula 

length 

metre 

m 

mass 

kilogram 

kg 

time 

second 

& 

electric  current 

ampere 

\ 

thermodynamic  temperature 

kelvin 

k 

amount  of  substance 

mole 

mol 

luminous  intensity 

c andela 

cd 

SUPPLEMENTARY  UNITS 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m/s 

activity  I of  a radioactive  source) 

disintegration  per  second 

(disintegration  |/s 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rads 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kg/m 

elec  trie,  capacitance 

farad 

K 

As  IV 

electrical  conduc  tance 

siemens 

S 

A/V 

electric  field  strength 

volt  per  metre 

V'm 

electric  inductance 

henry 

II 

Vs/A 

electric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

V/A 

electromotive  forte 

volt 

V 

W.A 

energy 

joule 

1 

N-m 

entropy 

joule  per  kelvin 

I'K 

force 

newton 

N 

kg-ms 

frequency 

hertz 

Hz 

(cycle)/s 

illuminam  e 

lux 

lx 

Im/m 

luminance 

c andela  per  square  metre 

cd/m 

luminous  flux 

lumen 

lm 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic  flux 

weber 

Wb 

V-s 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  fo  ce 

ampere 

A 

power 

watt 

VV 

)/s 

pressure 

pascal 

Pa 

N m 

quantity  of  elec  trn  it v 

c.oulomb 

<: 

As 

quantity  of  heat 

joule 

1 

N m 

radiant  intensity 

watt  per  steradian 

W'sr 

spec  if  if  heat 

joule  per  kilogram-kelvin 

I'kgk 

stress 

pascal 

Pa 

N/m 

thermal  conductiv  tv 

watt  per  metre1  kelvin 

W/m  k 

velocity 

metre  per  second 

m/s 

viscosity,  dynamic 

pascal  second 

Pa-s 

viscosity,  kinemati 

square  metre  per  second 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

| wavt’I'm 

work 

joule 

1 

N'm 

SI  PREFIXES: 

Multiple 

ation  Far  tors 

Prefix 

SI  Symbol 

1000  000  000  000  10” 

tera 

T 

1 000  000  000  10,# 

C 

1000  000  10* 

mega 

M 

1 000  io' 

kilo 

k 

100  10' 

hecto* 

h 

10  10' 

deka* 

da 

01  10" • 

deed* 

<1 

0 01  10- 1 

c «nti“ 

i: 

0 001  10  » 

mill* 

m 

0 000  (Mil  10  * 

micro 

u 

ooononnooi  io  v 

nano 

II 

0 OOO  OOO  OOO  001  10  ” 

pico 

p 

0 000  000  000  (100  001  10  •' 

femto 

f 

0 000  (HHMKH)  OOO  OOO  001  10  »■ 

alto 

a 

' To  be  avoided  whom  possible 
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